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METHOD FOR CALCULATING THE AERODYNAMIC LOADING ON AN 
OSCILLATING "FINITE WING IN SUBSONIC AND SONIC FLOW 
By Harry L. Runyan and Donald S. Woolston 


StMlAKY 


This paper presents a method for dete rmining the air forces on an 
oscillating finite wing of general plan form in subsonic flow including 
the limiting case of sonic flow. The method utilizes some of the concepts 
developed by Falkner (British R. & M. No. 1910) for steady lifting-surface 
theory. The loading on the wing is assumed to be given by a series con- 
taining unknown coefficients which satisfies various boundary conditions 
at the edges. The required integrations are performed by approximate 
means, and a set of simultaneous equations in terms of the coefficients 
in the loading series is obtained. Solution of this set of equations 
then gives the loading coefficients. The method is applied to rectangular 
and delta wings and comparison is made with existing theory. A sample 
calculation is given in an appendix. 


INTRODUCTION 


The analytical determination of the air forces on oscillating wings 
has been a continuing problem for over 30 years. Most of the effort has 
been directed toward the determination of the forces on wings in two- 
dimensional flow, and results have been obtained for a complete range of 
Mach numbers, both subsonic and supersonic. For finite wings, however, 
the analytical work is still in a state of development for all speed 
ranges. The main effort has been directed toward the incompressible 
case (for example, refs. 1 to 8) and the supersonic field (refs. 9 to 14) 
leaving an important area virtually untouched, namely, the subsonic range 
of Mach numbers between 0 and 1. Among the few studies that have been 
made in this range of Mach numbers have been those of Merbt and Landahl 
(ref. 15) , who extended some of the work at very low aspect ratio to 
include the effect of Mach number, and the work of Voss, Zartarian, and 
Hsu (ref. 16). In addition, Reissner (ref. 17) made a preliminary 
assessment of the use of Mathieu functions for subsonic compressible flow, 
and W. P. Jones (ref. 18) and Lehrian (ref. 19) studied the cases for 
larger aspect ratios. The present paper deals with some recent efforts 
which have been made toward filling this gap, namely the development of 
a wing-surface theory for subsonic and sonic flow. 
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In the stationary two-dimensional case, the effects of compressi- 
bility may readily be obtained from results of incompressible theory by 
application of simple transformation or correction factors such as the 
well-known Erandtl-Glauert factor. In the stationary three-dimensional 
case, the effect of compressibility may be treated as an extension of 
the two-dimensional procedure in a simple mann er as proposed by Gothert 
(ref. 20). In a further extension for the oscillating compressible case. 
Miles (ref. 21) has proposed a method involving several transformation 
steps for correcting the results of incompressible theory to include the 
effects of compressibility for a finite wing. The method is, however, 
restricted to low values of the reduced frequency. It therefore appears 
to be necessary to deal with the boundary-value problem directly in order 
to provide a method of general applicability which may be used, for 
example, in flutter work. 

The linearized boundary-value problem for the oscillating finite 
wing in a compressible flow may be approached from tiro points of view. 

One approach involves the transformation of the governing differential 
equation by a suitable choice of coord in ates followed by the use of the 
classical method of separation of variables. Solutions are then found 
in terms of series of orthogonal functions . This method was used by 
Schade and Krienes (ref. 22) in obtaining solutions for the oscillating 
circular plate and has recently been generalized by Kuessner (ref. 23). 

It has been pointed out that there is a definite limit to the number of 
generalized coordinate systems appropriate to the wave equation and to 
the corresponding orthogonal functions available. The method is, there- 
fore, limited to specialized plan forms such as the circular or the 
elliptic plate. 

The second approach involves a direct consideration of the integral 
equation relating pressure and downwash. The integral equation may be 
derived from the standpoint of either the velocity potential or the 
acceleration potential. Inasmuch as one is primarily interested in the 
pressures on the wing surface, it seems more direct to use the form of 
the integral equation associated with the pressure, namely, the acceler- 
ation potential approach; therefore, this procedure is used in the 
present paper. 

Basic to the solution of the integral equation is the evaluation 
of its kernel, a function which represents the downwash at a point in 
the plane of the wing due to a unit loading. Only recently (ref. 2k), 
the kernel function for oscillating finite wings in subsonic and sonic 
flow has been reduced to a form which can be conveniently evaluated, 

This background of available information on the kernel has led to the 
consideration of numerical approaches for handling the integral equation 
in the development of general lifting-surface methods which could be 
used for any speed range. 
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A similar integral equation appears in the steady case, hut 
involves a much simpler kernel function and a number of approximate 
solutions have been successful l y developed. In selecting an approxi- 
mate method for handling the integral equation in the osc ill atory case, 
the various procedures developed for the steady case were examined. In 
considering these procedures, it was kept in min d that a lifting- surface 
theory which is to be used for flutter calculations must meet certain 
requirements beyond those usua lly considered in the steady case. One 
essential requirement is that not only the lift but also the moment must 
be accurately predicted. In several of the methods now being used for 
the steady case, the wing is essentially replaced by a single line at 
the l/4— chord position. This placement essentially fixes the center of 
pressure at the l/4 chord, and thus the moment is constrained to become 
a function only of the lift force. In the case of an airplane of normal 
configuration with a tail placed a considerable distance from the wing, 
the contribution of the wing pitching moment is very small compared with 
the pitching moment due to the tail. However, for tailless configurations 
and for aeroelastic problems such as flutter, the moment on the wing is 
of predominant importance and must be accurately predicted. 

In addition to the previously mentioned requirement of predicting 
the moment, the method should be easily adaptable to the calculation of 
the loading on oscillating flexible wings and, in addition, should take 
into account the effects of compressibility. It appeared that, among 
others, a method developed for the steady incompressible case by Falkner 
(ref. 25) could be extended to include these requirements and this method 
was selected as the basis for the present investigation. 

Among other methods considered was a multiple-line approach suggested 
by the procedure of Schlichting and Kahlert (ref. 2 6) for the case of 
steady flow. The procedure is thought to be of interest in that it indi- 
cates a simpler alternative means of handling the lifting-surface problems. 
An extension of this method to the oscillatory case is also made in the 
present investigation. 

The primary aim of the present investigation has been the develop- 
ment of a lifting-surface method for calculating the forces on a wing of 
any plan form which is harmonically oscillating in a subsonic or sonic 
flow. The method uses some of the concepts of Falkner in an extension 
to the compressible oscillating case. First, a brief description of the 
basic integral equation and some of the possible means of solution are 
given. Then, the basic theory of the method is given in detail in the 
body of the report, and the results of some calculations are shown for 
rectangular and delta wings. In appendix A, the numerical details of 
calculating the forces on a rectangular wing are given as an example. 

In appendix B, a treatment is given for certain integrals which arise 
in the method and which contain singularities. In appendix C the expres- 
sion for the pitching moment for a delta wing is obtained. Finally, in 
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appendix D, a description and some results of applying a multiple-line 
approach, based on a procedure used by Schlicting and Kahlert for steady 
flaw, are given. 


SYMBOLS 


A 

a 


8 -rrni 





c 


F(M) 


a t 


f „(5> 

G < M > - m 


G ,G 
vo> vn 




aspect ratio 

axis of rotation measured from midchord in terms 
of half -chord, positive aft 

coefficients in expression for loading 

half -chord, ft 

lift coefficient associated with a degree of 
freedom q 

moment coefficient associated with a degree of 
freedom q 

chord, ft 


nondimens ional downwash factor defined by equation ( 29 ) 
chordvise loading function 


chordwise replacement loads 
spanwise loading functions 


K[M,|(x-|),|(y-Ti)] 

K(M,Z) 


kernel of three-dimensional Integral equation 


kernel of two-dimensional integral equation. 


K(M,Z) 



t> 


J 
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K nondimens ional three-dimensional kernel defined by 

equation (29) 

k reduced-frequency parameter based on half -chord, bo>/v 

4 

k reduced-frequency parameter based on e, k = 

L. total lift on wing, associated with a degree of freedom q, 

lb 

l reference length 

M Mach number 

M q moment about a, associated with a degree of freedom q, 

ft-lb 


n,m, j integers 


q 

R 1 = \[(x 
R 2 = \/U 

S 


8 


t 

V 

v(x,y,t) 

x,y,z,|,ri 


x e 


dynamic pressure, ipV 2 

- i ) 2 + P 2 (y + I ) 2 

- O 2 + e 2 (y - I ) 2 

area of wing surface 
semispan of wing, ft 
time 

stream velocity 

vertical induced velocity or downwash 
rectangular coordinates 
control-point location 


X Q = (x - Xy) 

x v 


y N 


chordwise coordinate for replacement loads 
distance from control point to center of segment 
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y Q - y - 


y 



R 


distances from a control point to midpoints of load segments 
located, respectively, to the left and right of the wing 
midspan, referred to e 


Z = k(x g - |) or k(x g - xy) 


a 


generalized coordinate of an angular displacement 



r(s,T|) 

AP v (p) 

Ap(s,n) 

6 

e 

9 

A 

P 

4 * 

0 


Subscript: 


vorticity distribution 

loading function defined by equation (26) 

pressure difference 
thickness ratio 

semispan of spanwise load segment, s/20, ft 
angular chordwise coordinate, x = -cos 0 
variable of integration 
air density 

angular spanwise coordinate, q = -cos 0 
phase angle between lift or moment and position 
circular frequency, radians/sec 


<1 


denotes a degree of freedom 


DISCUSSION OF INTEGRAL EQUATION AND METHODS OF SOLUTION 


The determination of the aerodynamic forces on an oscillating wing 
is a problem which, in general, must be handled by approximate or iter- 
ative procedures. In the two-dimensional case one approach to the 
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problem has been the transformation of the governing differentia] equa- 
tion and its associated boundary conditions to an integral equation for 
which approximate solutions can be found. It is interest ing to note 
that for this two-dimensional case, an iterative solution based on the 
integral-equation approach preceded a more exact solution based on a 
series of Mathieu functions by about ten years. 

In the case of wings of finite span, a transformation similar to 
that made in the two-dimensional case and leading to an ±ntegr|il equa- 
tion can be made. The integral equation may be derived on the basis 
of the velocity potential or oh the basis of the acceleration potential. 
The advantage of the use of the acceleration potential rests in the fact 
that the trailing surface of discontinuity is not explicitly introduced 
into the problem and the integral equation relates directly the downwash 
and the pressure difference on the wing. The disadvantage of the accel- 
eration potential as compared with the velocity -potential approach lies 
in the appearance of a more complicated kernel of the integral equation. 


Integral Equation 

For compressible flow, the integral equation as derived on the 
basis of the acceleration potential and given, for example, in refer- 
ence 2k may be written as 


w(x,y,t) 

V 


1 /TApCl^T]) 

EjtJJ q 


K 


M,|(x-£),^(y-q) 


d£ dq 


( 1 ) 


where w is the known vertical induced velocity or downwash, Ap(|,q) 
is the loading or pressure difference (positive upward) at a point (|,q) 
on the wing surface, V is the stream velocity, and K is the kernel 
of the integral equation. The pressure difference Ap(|,q) is related 
to the bound vorticity through the equation 


Ap(| ,q) = pVT(£,q) 


so that equation (l) is sometimes written with 

2T(|,q) 


Ap(|,q) 


replaced by 


V 


The function K represents, physically, the contribution to 

.the downwash at a field point (x,y) in the plane of the wing due to a 
unit pressure difference at a loading point (£,q) also located in the 
plane of the wing. 
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The kernel K is discussed in reference 2k and appears originally 
in an integral form as 

i-^ A-M\/ A 2 +p 2 (y-T) )2+p 2 z 2 

^ e -^)' 


KM,|(x-|), f §(y-T\ ) = 11m 

L v v J z — > 0 


Vp* 




<5z c 


s yA 2 +p 2 (y-Ti) 

where to is the circular frequency of the oscillating wing 


<37v 


2 +P 2 z 2 


( 2 ) 


In the form given by equation {2), the kernel function appears as 
an improper integral whose complete Integration in closed form has not 
been accomplished. The kernel function has been extensively treated in 
reference 2k, however, and has been reduced to a form which can be eval- 
uated. The reduced form of the kernel given in reference 2k contains 
Bessel and Struve functions and proper nonsingular integrals which can 
be handled by numerical means. The kernel contains singularities at 
y - 11 = 0 and x - £ > 0, but these singularities have been isolated 
and expressed in a form which can be handled in numerical procedures. 

Also given in reference 2k is an expansion of the kernel into a series 
in powers of the reduced frequency, which for the present time serves 
as the most practical form for application of the kernel. 

The solution of equation (l) requires a determination of the unknown 
loading distribution Ap(£,n) subject to known boundary conditions. In 
view of the complicated nature of the kernel, an exact analytical solution 
of equation (l) does not seem possible. It appears necessary, therefore, 
to consider some numerical method of solution and several possible approx- 
imate methods are discussed in i)he next section. 


Tf 


Some Approximate Procedures 

In searching for an approximate means of handling the integral 
equation, two approaches seem likely, one involving an approximation 
for the kernel and the other an approximation for the loading. Hie first 
of these would involve a replacement of the complicated kernel function 
by a simpler function which behaves In a similar manner (with regard to 
singularities, and so forth) but with which the Integral equation could 
be Inverted and solved analytically for the unknown loading. Such a 
procedure has been used by Fettis (ref. 27) for the two-dimensional, 
subsonic, oscillatory case and, in part, by Lawrence (ref. 28) for the 
finite wings in incompressible steady flow and by Lawrence and Gerber 
(ref. 29) for oscillating finite wings In incompressible flow. The 
application of such a procedure to the case of compressible flow would 
be quite difficult, but is perhaps worthy of further consideration. 
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A second possible approach involves leaving the kernel unchanged 
and selecting an appropriate expression for the loading which leaves 
only certain unknown coefficients to be determined. By such a procedure 
the integral equation can be reduced to a sum of def init e integrals anfl 
the unknown constants can be determined by collocation, that is, 'by- 
forcing the solution to fit the known downwash conditions at a number 
of points on the surface. Although difficulties arise in the integra- 
tions because of singularities in the kernel and the loading, they can 
be handled in numerical methods. Various forms of this procedure have 
been used quite successfully for oscillating two-dimensional wings anri 
for finite wings in steady flow. 

A logical extension of this second approach, the choice of a 
replacement function for the loading, to the three-dimensional problem 
would involve finding a means of handling the definite integrals which 
arise. Xt would seem possible to divide the wing into many smal 1 areas 
and, with an assumed form of the loading on each area, integrate numer- • 
ically over each area to find the downwash at a point on the surface. 

Such a procedure would lead, however, to a system of equations in the 
unkn own constants of at least the order of the number of areas so that 
lengthy computations would arise which would require computing equipment 
of great capacity. 

In searching for a less cumbersome method, it is natural to study 
the possibility of utilizing some of the concepts of procedures which 
have already been developed for the case of steady flow. Several methods 
for steady flow seem adaptable to the unsteady problem. Among those to 
be considered are the procedures of Falkner (ref. 25), Multhopp (ref. 30), 
Weissinger (ref. 31) , and Schlichting and Kahlert (ref. 26). A brief 
discussion of these methods in relation to the objectives of the present 
investigation is given subsequently. 

The methods developed for the steady case by Falkner and by Multhopp 
represent attempts at lifting-surface procedures for handling the integral 
equation. In both approaches, the unknown loading is expressed in terms 
of a series of chosen modes of loading with unknown coefficients to be 
determined by collocation. In both methods, also, the double integrations 
are performed numerically and lead to sets of downwash factors. It 
appears that either method could be extended to the unsteady case nnfl 
systematized through the development of suitable tables of the downwash 
factors. For the purpose of extending such procedures to the unsteady- 
flow problem, in the present investigation, the Falkner procedure has 
been considered preferable since considerably fewer tabulations are 
required. The development of a lifting-surface approach based on the 
method of Falkner is the main objective of this paper. An extension 
of the steady-state procedure of Multhopp to the oscillatory case has 
been performed by Jordan and some results are given in reference 32. 
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It should be noted that W. P. Jones (refs. 2 and 18) has also 
treated the problem of determining the aerodynamic forces on wings of 
any plan form by usin g some of Falkner * s concepts. The method differs 
from the one to be presented in this paper and pursues a different numer- 
ical pa th since in the Jones procedure the basic integral equation was 
treated on the basis of the concept of the velocity potential. The use 
of the velocity potential results in a double integration over the entire 
field of velocity discontinuity, that is, over the wing and its wake. 

The use of the acceleration potential, employed in this paper, involves 
a double integration only over the area of pressure discontinuity, that 
is, over the wing surface itself, since no pressure discontinuity exists 
in the wake . 

Other steady-state procedures mentioned previously were those of 
Schlichting and Kahlert and of Weis singer . These methods represent 
lifting-line rather than lifting-surface approaches. In the Weissinger 
method a single lifting line is placed at the l/4- -chord position on the 
surface and the known downwash is satisfied at points on the 3 A -chord 
position. Th e Schlichting procedure makes use of several lifting lines. 
The surface is divided into a number of spanwise strips of equal chord. 

A lifting line is placed at the l/4- chord position of each strip and the 
downwash conditions are satisfied at points on the 3 A -chord position 
of each strip. Since a lifting-line procedure might be expected to 
involve less labor than a lifting-surface approach, an attempt was made 
in the present study to extend the Schlichting method to the case of 
unsteady flow. The application of the Schlichting approach to the 
oscillating case is discussed in appendix D. 

In the following two sections the lifting-surface method of the 
present paper is developed, first for the case of subsonic flow and 
then for the limiting case of M = 1. 


DESCRIPTION OF SURFACE -LOAIENG METHOD FOR SUBSONIC FLOW 


The solution of the integral equation for the case of an oscillating 
three-dimensional wing involves a double integration over the plan form 
of the wing as indicated by equation (l): 


w(x,y,t) 

V 


1 

6kq 



Ob \ <b. . 

M,-(x-|),-(y-q) 


d£ dt) 


Of the three main ingredients of the problem - the downwash, the kernel, 
and the loading - which are related by this equation, two can be considered 
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a6 known. The kernel function K can be evaluated through the use of 
the forms given in reference 24. The downwash w(x,y,t) is determined 
by the motion of the airfoil. This leaves as the only unknown in equa- 
tion (l) the function Ap(|,q), the magnitude of the loading at points 
(t,T|) of the surface. 


Statement of the Boundary Conditions for the Downwash 

For linearized flow, the downwash function v(x,y,t) can be deter- 
mined from the motion of the mean surface of the airfoil thr ough the 
relation 


w(x,y,t) 


Jj: f ( x >y^) + V ^ f (x,y,t) 


(3a) 


where, for harmonic motion, 


f (x,y,t) = f (x^e 10 * 


so that f (x,y) is the amplitude of the displacement of the mean surface 
from the equilibrium position. Ctnce an appropriate function for f (x,y) 
is chosen, the downwash is known and the loading can be dete rmin ed from 
equation (l). For the usual modal type of flutter analysis, for example, 
f (x,y) could be found from the assumed vibration modes of the structure. 

Because the basic problem of determining the loading for a given 
downwash as expressed by equation (l) is linear, it is often convenient 
to consider the problem separated into a sum of individual problems, 
each associated with a particular selected type of motion or degree of 
freedom and, hence, with a particular downwash function. This implies 
that the total downwash can be written as a sum of individual downwash 
functions, so that equation (3a) appears as 


w(x,y,t) 



f x ( x ^y) + f 2 ( x >y) 


+ 


+ 


f q ( x ,y) 


hut 

e 




(3b) 


where w^(x,y) is the downwash associated with the particular degree of 
freedom f^(x,y). The basic problem of this paper is then the determina- 
tion of the loading Ap^(|,T|) associated with any one of the downwash 
factors Wq(x,y) . 
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Treatment of the Load in g 

La addition to the downwash conditions set "by equations (3), the 
loading Zg?(|,T]) must satisfy various conditions at the edges of the 
plan form. In the chordwise direction, for subsonic flow, the loading 
must he infinite at the leading edge and vanish, in accordance with the 
Kutta condition, at the trailing edge 0 In the spanwise direction, the 
loading should become zero at the wing tips. 

Series expression for the loading .- As a first step in the deter- 
mination of the unknown load in g, it is assumed that the function Zp(|,T]) 
can be expressed in terms of a series of both spanwise and chordwise 
pressure or loading modes, so chosen as to satisfy the edge conditions 
just discussed and containing arbitrary coefficients to be determined. 

The spanwise modes are expressed in terms of the variable T|. It is 
convenient to express the chordwise pressure modes in terms of an often 
used variable 0, related to £ and based on the circular cylinder 
enveloping each chord as diameter and shown by the following sketch and 
in the following equation: 



6 = fmlOl) “ h(q)cos 0 (4) 

where 

b 0l) = l^te^) “ f *e(n)] 

In terms of the variables 0 and T], the loading associated with 
a particular downwash distribution w^(x,y) may be expressed in terms 

of the s ummat ion of pressure modes and written as 
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v s ' ,|) = 


4 


, 2 2 
S c - Tp 


If 


cot 2-(a^ + a^Ti + a^ D ri 2 + 


00 


02 



sin 0 



q. 

V + 



+ sin 20 



q. 

+ a 2 1 'n + 


<1 

a 22 n 



(5a) 


where the superscript q identifies a particular downwash. 

The series form of the loading given by equation (5a) satisfies 
the Kutta condition (Ap = 0 at the trailing edge) and has the desired 

type of singularity 1 / sjx wl lere x is measured from the leading edge. 
The terms in T) are of such a form as to cause Ap(|,i]) to become zero 
with infinite slope at the wing tips. The use of this series form for 
^>(i,n) means the pressure distribution is essentially synthesized by 
a series of chosen pressure modes and that the values of determine 

the contribution of each pressure mode to the final pressure distribution. 

Equation (5a) can be written in more concise farm, as 


V 9 ’ 1 ') = 1ST ^ (5b) 

11=0 m=0 

where 

f 0 (O = cot | 
f n (|) = sin n0 
^(q) = q m \/s 2 - q 2 

(For symmetrical motions of a wing, w <1 (x,y) =1 w^(x,-y), the odd 

terms in q need not be retained. For unsymmetrical motions’, 

Vq(x,y) = -Wq(x, -y), the even-powered terms in q do not have to be 

retained. In the numerical examples to be given in a later section, 
symmetric motions are considered. ) 


(n = 0) 
(n / 0) 


J 


(6) 
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Illustrative expressions for the force and moment .-* The specific 
problem is the dete rminat ion of the coefficients % in equation (5a) 
or (5b) and, as has been indicated, a particular set of values for a m 
must be found for each downwash v^(x,y) . Once a set of values for a^ 
has been determined, a particular loading Apq(|,Ti) is defined by equa- 
tion (5a) or (5b) • This expression for Zp^(|,Tj) gives the pressure 

distribution and can be integrated to obtain various desired section 
forces or total forces. For_example, total lift and moment coefficients 
associated with the motion f^(x,y) can be found through the relations 



(7a) 


and 



1 

1 ~~ qSb 



a)Ap q (!,n)<3.£ < 1 ti 


(7b) 


•where the moment is obtained about a spanwise axis | = a. 

By substituting equation (5a) into equation (7a) the following 
result for the lift coefficient is obtained which is valid for any 
plan form: 



* 2 A [ r <1 1 .1 1 1 O 

(l6 a oo + ° a io + AaQg + 2a;i_2 + 2a o4 + a l4. 


(8) 


The expression for the moment about the axis considered must be 
determined for each particular plan form. The moment coefficient will 
take on a form similar to equation (8) but with constants determined 
by the plan-form geometry. For a rectangular wing, for example, the 
moment coefficient for moment about the midchord (a = 0) is 



A 

52 


16a 


00 


1 1 1 

+ ka.^ + 2 aQ4 


q. 

*20 


+ 8a + 2a 


q. 

i 22 


(9) 


The equation for the moment coefficient for a delta-wing plan form 
is given in appendix C. 
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Replacement of the Integral Equation by a S umma tion 
of Definite Integrals 


In the rest of the analysis, it will be convenient to deal with a 
general downvash w(x,y) and to drop the subscript and superscript q. 
In order to evaluate the coefficients a^ for a particular downwash 

w^(x,y), the first step involves the substitution of the loading series 
(eq. (5a)) into the integral equation (eq. (l)) to give 


w(x,y) _ 4pV 2 
V fircqb 


4 


2 2 
s^ - T) C 


cot |(a 00 + a 01 n + a 02 T ) 2 + ...) + 
sin e(& 1Q + ^T] + a^q 2 + . . + sin 26 (b.^ q + a^i) + 


a 22^ + 




— — 

) + • . . 

K 

M,^(x-|),|(y-ii) 


?L6 dq 


( 10 ) 


or if the form for Ap(|,T|) given by equation (5b) is used, the following 
is true: 


w 


^ II tl v JJ^t )%,(i )] « ^ (id 

11=0 m=0 a L J 


A significant step toward the desired solution of the integral 
equation has now been arrived at because equation (10) or (ll) Is no 
longer an integral equation, that is, an unknown quantity no longer 
appears under an integral sign. Instead, the downwash w(x,y) is now 
given by a summation of definite integrals multiplied by various unkn own 
constants a^. These integrals are of the form 


JJ cot | r\ m sjs 2 - T ] 2 k[m,^(x-| ),®(y-ri)J d| dr] 

S 

JJ sin nG r\ m \js 2 - n 2 K[M,|(x-|),^(y-Ti^j d| ^ 

S L - 1 




( 12 ) 
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The numerical evaluation of these integrals is an important part 
of the labor of the present method, for once their values have been 
calculated, one can proceed immediately to the determination of a^ 

and, therefore, to the determination of the pressure distribution. In 
the method under consideration, the required integrations are performed 
numerically. Because the chordwise and spanwise integrations are 
treated differently, they are discussed separately in the next two 
sections, the chordwise integration being discussed first. 


Treatment of the Integration in the Chordwise Direction 

For the purpose of performing the chordwise integration it is 
convenient to replace the right side of the integral equation by a sum 
of separate integrals, each involving a particular chordwise pressure 
mode and in which the spanwise (rj) and chordwise (| ) integrations are 
separated. For this purpose equation (ll) may be written as 


w(x,y) 


1 kpv^y*- 

b m=o 


I 


a 0m J "** 'l a nm / 

span n=l span 


( 13 ) 


where 


Jo= f f 0 (S )k|m,~:(x-£ ) j®(y-n )1 d| 

J chord L v v J 

and 



(14-a) 


(14-b) 


Remarks on methods of performing the chordwise integrations. - One 
could at this point evaluate the integrands of Iq and ^ at small 

intervals of g and integrate by numerical means. This procedure would 
be difficult and tedious, however, since the kernel contains singularities 
and involves a large number of parameters and since the integrations would 
have to be carried out for many values of t|. It is desirable to seek a 
simpler, but more approximate, approach which would be less laborious. 
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3h the present analysis, an approach developed by Falkner in ref- 
erence 25 (and discussed more f ully in ref. 33) for the steady case is 
used -which greatly reduces the amount of labor. In essence it is assumed 

that in the integrals Iq and Iq, the loading modes fQ(|) = co "^ 2 ar ^ L 

f n (| ) = sin n0 can be replaced by an arbitrary but small number of 

replacement loads which are chosen in order to satisfy some integrated 
properties of the continuous loading. The use of this assumption w ill 
result in the integrals Iq and Iq appearing finally as summations 

In the forms given by equations (20a) and (20b). 

The use of a set of Individual loads in performing the chordwise 
integrations and the procedures for evaluating them introduce approxi- 
mations which are somewhat arbitrary. These approximations anri their 
implications will be discussed later as they arise. They correspond 
to the assumptions made by Falkner and have been used in the present 
investigation for practical reasons. With high-speed computing equipment 
available, some of the approximations in the development which follows 
might be avoided. It can be stated, however, that, in the cases consid- 
ered, the end results seem to be satisfactory. 

Calculation of the replacement loads .- In calculating the values of 
the replacement loads, it is first required that the sum of the individual 
loads associated with a particular chordwise pressure mode must equal the 
integrated value of the continuous chordwise loading. If the continuous 
distribution is replaced by j replacement loads and G^q and Gv n 
denote loads associated, respectively, with the modes cot 0/2 nnrl 
sin n0 , this equivalent load condition leads to the following set of 
equations : 


o / 0 

Z_ G v o = / co ' t o 6111 9 d0 = it 
v=l J 0 


(15a) 


and 


T~ G = / sin n0 sin 0 d| = ^ (n = l) 
v=l v JO 



> (15b) 


V=1 


(n > 1) 
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The use of equations (15) imposes one condition on the replacement 
loads G. However, equations (15) contain j unknowns; therefore, at 
least j - 1 additional equations in terms of G are needed. In the 
present method the additional equations are obtained by requiring the 
downwash produced at selected points on the chord by the individual loads 
to equal the downwash due to the continuous load distribution. In 
imposing this condition the two-dimensional kernel function K is used 
Instead of the more complicated three-dimensional downwash factor defined 
by equation (2) . The relation between the two- and three-dimensional 
kernel functions is given by the following equation: 


k[m,®(x-i)] =J K|M,ffi(x-6),“(y-n)jdTi 

where K [kf (*-!)] is the two-dimensional kernel. 

The use of the two-dimensional downwash factor at this point is 
an intermediate step used only for evaluating G. In satisfying the 
actual downwash on the wing in the next section the calculated replace- 
ment loads are used with the three-dimensional kernel function. This 
intermediate use of the two-dimensional downwash factor implies that 
the areas close to a point contribute the major portion of the downwash 
at the point. 

The j - 1 additional equations in terms of G can be written by 
selecting j - 1 positions, or control points, on the wing chord denoted 
by Xg at which to equate the downwash produced by the individual loads 
to that produced by the continuous loading. Although the lift conditions 
(eqs. (15)) did not require any assumption regarding the location of the 
replacement loads, for the purpose of employing the downwash condition 
such an assumption is necessary. In selecting this location an inter- 
pretation of the meaning of the replacement loads G is somewhat arbi- 
trary. The replacement' loads G can be considered as loads distributed 
over a certain portion of the wing chord, or they may be considered as 
concentrated loads operating at particular points on the chord. If the 
concept of a load distributed over an area is taken, then the distance 
to the downwash point is physically identifiable as some average distance, 
which for this particular case is measured from the center of the load 
area. If the concept of the concentrated load is used, then the distance 
is measured from the point of application of the load. In either case, 
a reference point x v is selected as the point of application of the 
load G so that the distance to the downwash point becomes Xg - x^. 

This distance is then used in developing the j - 1 downwash equations 
for the determination of G. 
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For the cot 8/2 term of the loading series, the downwash equations 


El = f oot a sin 6 

V=1 o-l ^ 


(l6a) 


and, for the loading sin n8, the downwash equations became 


_L _ ol 

> GynKgv = / Bin n0 slri 9 ^gl de 

v=l J -1 


where one such equation is obtained for each control-point location 
Xg selected. The functions Kgy and Kg^ in equations (16) are 

defined as follows: 


%; = K(M,Z„) 

^ and 

- E ( M ’ z i) 

and represent the two-dimensional downwash factor tabulated by Schwarz 
(ref. 3*0- I 11 these functions Zy = k^xg - x v ) and Z| = k^x g - 

where k = — and where a semichord b has been chosen for convenience 
V 

as a reference length. 

It should be noted that the integrands in equations (l6a) and (l6b) 
become infinite at £ = x and the integrals must be carefully treated. 

The method of handling these integrals is given in appendix B. 

Other conditions which would not bring in consideration of the two- 
dimensional kernel could presumably be used for determining values of G. 
For example, in place of using the downwash conditions, j - 1 additional 
equations could be written for each chordwise term of the loading series 
in which the first and higher (through j - l) moments of each load G 
could be equated to the corresponding moments produced by the continuous 
loading. By such a process, however, the loads G would not be a func- 
tion of either frequency of oscillation or of Mach number. Through the 
“ use of the downwash conditions, some effects of both frequency and Mach 

number are included at this intermediate stage. 
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'niustrative equations for the loading functions . - In applying the 
method of this paper, three chordwise pressure modes have been retained. 
The continuous distribution in each mode has been replaced by four 
Individual replacement loads. 

In order to evaluate the replacement loads, the lead established 
in the steady case -with regard to locations of control points and loads 
has been followed, Hie loads are assumed to act at the l/8-, 3/8-, 5/8-, 
and 7/8 -chord positions, and the downwash conditions are applied at 
control points midway between the load stations, that is at the ifk--, 
l/2-, and 3/k—chor& positions. Replacement loads based on these chord- 
wise positions have been calculated for M = 0, 0.5* and 1.0 and are 
presented In table I. 

One set of replacement loads must be determined for each chordwise 
pressure mode. For the mode involving cot 0/2, the values of 

are determined from the following equations: 

G 10 + G 20 + G 30 + So = f 0 

G loSl + G 2 oS_ 2 + G 3oS.3 + 

SoSl + G 20^22 + G 3(fe + 


G 10 K 31 + G 20 K 32 + g 30 k 33 + 


cot g sin. 0 d0 = it 


^ So 


cot sin 0 d0 


> U?) 


SO -ot|sin0K 2| 


d0 


^ 0 % =f 


cot 2 sin 6 d0 


where the functions K used here and In equations (l8) and (19) are 
defined after equation (l6b). 

The first of these equations imposes the lift conditions stated by 
equation (15a) . The remaining equations apply the downwash condition 
of equation (l6a) at each of three control points; the second equation^ 
for example, states that at control point 1 (at Xg = x^ in eq. (l6a)J 

the downwash produced by the four loads at Xy_-^, x v _2* x v~3* aQ ^' 
Xj,-^ must equal the downwash produced at control point 1 by the con- 
tinuous cot 0/2 loading. 
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A similar set of equations can be -written for each chordwise pres- 
sure mode by using equations (l^b) and (l6b). For sin 8 the set of 
equations is 


G n + G 2 i + G 3i + ~ f 0 sln2e de " I 

°11%1 + g 21%2 + ^l^OS + %l K lA- = f 8 11129 K l| d0 

JO 


G 21 K 21 + g 21 K 22 + ^1^23 + ^l 1 ^ 
+ G 21 K^ 2 + G 3l^3 + fyl^ 
and for sin 20, 



sin^Q d8 
sin 2 0 K^| d0 


■\ 


G-] p + Gpp + G^2 t 

G 12 K I1 + G 22 K 12 + 
G I2 K 21 + G 22 K 22 + 
G 12%1 + G 22 K 32 + 


-=r— 


G 32 K 13 + %2 K l4- 

G 32^23 + ^ 2^4 
G32&33 + Glj .^%4 



sin 0 sin 20 K-j^ d0 

> 

sin 8 sin 20 d0 
sin 0 sin 20 K^g d0 


(18) 


(19) 


Use of the replacement loads G in the chord-wise integration .- Once 
the replacement loads G have been determined, by solving equations ( 17 ), 
(l8), and ( 19 ) , continuous integrations indicated by equations (l4) can 
be replaced by summations of the products’ of the loads G and the three- 
dimensional kernel function as follows: 
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^> = b 


T~ G v0 K ^^ x " x v)^f(y- T i)l 

V=1 L- - 1 


and 3^ (eq. (l4b)) becomes 


r n = *TZ. G vn K t^( x - x v)^(y-n) 
v=l L 


(20a) 


(20b) 


"where the factor b arises through use of the substitution £ = -b cos 0. 
In effect the continuous integrations have been represented by the 
replacement loads and a mean value of the kernel. Ihese expressions 
may be substituted into equation (13) as follows: 


w(x,y) VV 2 sr~ r 

2_ 2 2 


v mS 


vn 


f gmCn)^ M ffi(x-x v ),®(y-Ti)ldTi (21) 
'■'span L v » J 


leaving only the spanwise integration to be performed. 


Spanwise Integration 

For the purpose of discussing the spanwise integration, equation (21) 
may be rewritten in simpler form as follows: 


>r( x >y) = _1_ 

V fkq 


n / 

V=1 


span 



( 22 ) 


where 


AP y (q) 


00 00 


4pV^ 'y a rrrn^yri6m (*1 ) 

n=0 m=0 


(25a) 
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or 


AP v (n) = \Jb 2 - n 2 GyO^O + W + * 02 ^) + ^vl ( a 10 + 

a ll T l + a 1 P 1 ! 2 ) + ^v2 (^20 + a 21 T l + a 22 T l^ + • • •! 


( 23 b) 


Equation (23b) may be seen to correspond to the original form of 
the loading Ap(i,T}) given by equation (5a) . In equations (23), how-, 
ever, the continuous chordirf.se -loading terms cot 0/2 and sin n9 have 
been replaced by the loads G v q and G yn . 

The integration of equation (22) may be handled by several proce- 
dures. A straightforward numerical integration could be performed for 
each value of x v by evaluating the kernel K at a number of spanwise 
stations. However, the kernel K contains some difficult singularities 
which have to be carefully considered. Moreover, in order to male tables 
of these integrals for general use, tables would have to be made for 
every aspect ratio, sweep angle, Mach number, and reduced frequency. 

In order to facilitate the development of tables, it is desirable to 
make use of a mean-value integration which involves integrating the 
kernel over a short segment of span and using the value of AP V at the 

midpoint of the segment. Th i s method was followed in the present anal- 
ysis and has the advantage that the integrals have to be tabulated only 
as functions of Mach number and reduced frequency. 

The downwash at a point (x,y) due to a small element of span of 
length 2e, the center of which is located at | = Xy and q = is 


where 


Aw(x,y) 

V 




(2k) 


(25a) 


and where AP^y^ refers to the value of AP y at the midpoint of 
the span element over which the integration extends. 
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The integral given by equation (25a) is of central importance in 
the present method and represents a "downwash factor" giving the down- 
wash at a point (x,y) due to a unit pressure loading acting over a span 
element of length 2e. The value of the integral dep end s only on the 
relative distance of the element from the point (x,y) and not on the 
spam-rise location of the element on the wing. For convenience of cal- 
culation, therefore, it is desirable to perform a coordinate transfor- 
mation to the center of the element, so that the integral appears as a 
function of distances from the point (x,y). 

In order to perform this transformation, let T| * = r\ - y K in equa- 
tion (25a) where rj * is a new spanwise variable. Equation (25a) c an 
then be written as 



k[m,|(x-x v ) ,f(y-y N -V)] av 

k |m,^(x-x^ ,^(y'-T]*)j dq* 


(25b) 


where y* = y - y^ and where y’ and V are shown in the following 
sketch: 
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Selection of Reference Lengths 

It is convenient at this_point to choose certain reference lengths 
In both the downwash factor F v jj, defined by equation (25b), and in the 

loading AP V (y-^ , defined by equations (23). In the case of the loading, 

the wing semispan s is chosen as a convenient reference length. If a 
new variable p = y^/s is introduced, the loading AP v (y N ) can be 

written as 


AP v (p) 





(26) 


where 


5n(n) = 


and where p is measured to the midpoint of the spanwise segment. 

In the case of the downwash factor F v ^, the length e is taken 
as a reference length. The variables x, x v , y 1 , and rj ’ are then 
considered in a new sense to mean nondimens ional quantities obtained by 
dividing_the dimensional distances by the length e. The downwash 
factor F vN then appears as 



where 


k = 


ecu 

V 


The use of e as a reference distance introduces a factor l/e^ 
in the kernel (eq. (2)) as follows: ' 
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K|il,k(x-|),k(y-Ti)J = 


_ ox-i 

J2_ -ttfcH) / £ 


^|[a-h x// 2 +e 2 (y-t! i 2 +e 2 i 2 ] 


11m “-5- e / 

z — ^ 0 dz J _ 




SK> ( 27 ) 


Equation (27) can be used to define a nondimensional kernel function K 
by writing 


K^k^x-x^ ,k(y’-V)j = K ji^k^x-x^ ^(y'-q 1 )] (28) 


so tbat £ in equation (28) corresponds to the quantity in braces in 
equation (27). 

With the use of equation (28), the downwash factor P v jj appears 
as 

F vN = | J' K^kpc-Xj,) , k(y 1 -q ' )j dr] ’ 

Finally, a nondimensional. downwash factor F v jf can be defined as 


F vN = 



(29) 


Final Expression for the Integral Equation 

Equations (2 6) and (29) can now be substituted into equation (24) to 
write a final expression for calculating the downwash at a point (x,y) 
on the wing. If the downwash due to all - the span elements over the 
wing is summed, the total downwash becomes 


w (x,y) 


1 

8rtqe 


£ X 

V=1 yw 


APy (M') F vW 


V 


(30a) 
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Ey using the expanded form for AP v (p) given by equation (2 6), the 
total downwash can be -written as (for e = s/20) 


w(x,y) 
— T~ 


n m v 


y s mn G vn®m ^ ) -^vN 

y N 


(30b) 


The only unknowns in this equation are the coefficients a^,. Ey 

selecting a number of control-point locations (x,y) at least equal to 
the number of a^'s desired, equation (30b) can be written for each 

control point. Thus, a set of simultaneous equations is obtained in 
terms of the coefficients which can then be calculated. 

The next section is concerned with the calculation of the downwash 
factors F vN which have been defined in equation (29). 


Determination of Downwash Factor F 

The value of the downwash factor F (eq. (29)) for the oscillatory, 
compressible case is not available in the literature and consequently 
must be derived. However, the kernel K has been discussed in detail 
in reference 22, in which several forms of K are given. One form, as 
given by equation (20) of reference 22, would require a numerical inte- 
gration in order to obtain F, since K cannot be integrated in closed 
form. Another form of the kernel, an expansion in terms of the frequency 
parameter k, is given by equation (5*0 of reference 22 and is the form 
employed in the present analysis. In reference 22, a discussion is given 
of the accuracy of the frequency expansion. It is indicated that for 
moderate values of E and M the series is quite good for small 
distances x - x v and y 1 - T] 1 . As Tc and the distances increase, 
the series deteriorates. Overall accuracy is not seriously affected, 
however, since the major contribution to the downwash at a point is 
from the nearby points on the surface. 

Accordingly, the series is integrated term by term to obtain an 
approximate value for F. Several of the integrals are, however, 
improper, due to the existence of s ing ularities in the integrand, and 
the concept of the principal value or finite part must be utilized. 

(See, for example, ref. 35 for a discussion of the finite part of 



Integrals.) The results of the Integration for an expansion to the 5th power of k is given in 
equation (51) . (in presenting the expression for the downwash factor F, it is convenient to drop 
the primes on the quantities x' and y’ and the bar on the reduced-frequency parameter k. In 
the following expression, therefore, x. y, and k will be considered to denote the dimensionless 
quantities referred to e unit length.) 


FyH - J~_ K(Wl) d 1 


-ita 0 / *0 + % *o_ + *2_ J* P (y - + *2 


— ■ — — — — ■ 

x Q (y + 1) x Q (y - 1) 0 P(y + 1) +'^1 2p 

+ P 2 (y + 1) 


+ -~[2P 2 + 34 - 2p 2 log + 


*''> - «h ' * % 


Xq + Rq 

r 6 FFT ij 


to? 


*o + 


log ?; — r r - log 
3M 2 x 0 2 


k(y + 1) 


1 


i - P (y - 1) .+ R 2 , 
2p OB p(y + 1) + Rq 


3M 2 - 2 


[(y - 1)R 2 - (y + l)Rq]| + ^,-33^x 0 2 - | MP 2 (5M 2 - 3)(3y 2 + l) - 

J 9&P L 


x Q + R 


20 6 (y - l) 3 log *° 7 -— + 2p 6 (y + l) 3 log 

(M + l)(y - 1) (m + l)(y + l) P 


h o \ 3 p(y - 1) + R 2 

— (m 2 + 4) x^qog — — - — + 

° v 1 p(y + q) + R]_ 


(2 - 1«4 2 - 4K 4 ) Xo \{j - i)Rg - (y + i)Rqj + P 6 [(y + I) 3 - (y - (-27 + | + | - *i) 

2 p 6 (y - l) 5 log - 20 6 (y + l) 5 log + M 2 )xq 3 - 


ro 

00 


(Equation continued on next page) 


*- 


•fl69C fCG VOVH 



l» A "I \ . Tl 

HU - -r ±\2 


+ ij)xo^log 

8p ^ p(y + l) + Ri 


- &pfi 2 (3y 2 + i)x o - ^[(y - 1) 3 R 2 - (y + 1)\J + 

|(6 - 2CM 2 - 15M 4 ) [(y - l)Rg - (y + l)R^ x Q 2 + ^(sm 2 - l) T(y - l)R 2 5 - (y + l)Ri^j - 
l8p®y(y 2 + l) + 2f3 8 [fr + l)^ - (y - 1)^ 


(31) 


where 


Ri = ^Jx 0 2 + P 2 (y + l) 1 


r 2 = \/ x o 2 + p 2 (y - 1) 2 


*o 


B. X - 


7 = 0.5772157 


where Xq = x - x v is the distance from the line of integration to a control point, made nondimen- 
sional by dividing by e, and is positive for locations of the control point behind the line of 
integration. The distance y is also nondimens localized by s and is positive to the right and 
is measured from the center of the integration segment as previously discussed. Ihe reduced 
frequency k must alBo be baBed on the span segment e and is accordingly 


k 


ecu 

V 


Of course the number of spanvrise segnents N can be arbitrarily chosen. In general, the more 
segments taken the more accurate the result. F alkn er (ref. 25) has used e <=< s/20 for most cases. 
As an example, the layout for a rectangular wing as used in the numerical example in appendix A is 
shown in figure 1. 

The next section of the paper is concerned with the application of the surface-loading method 
to the special case M - 1.0. 


ro 

vo 
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THE SURFACE-LOADING METHOD FOR SONIC FLOW 


The problem of calculating the forces for the limiting case of a 
wing -which is oscillating in sonic flow is in general the same as for 
a wing at subsonic speeds with two exceptions. First, for the case of 
a sonic trailing edge, that is, when the trailing edge is perpendicular 
to the flow, the satisfaction of the Kutta condition (Ap = 0 at trailing 
edge) is no longer necessary and another form of the series expansion 
for the loading Zip is applicable. (For the subsonic -trailing-edge 
case, i.e., when the tra iling edge is not perpendicular to the flow, the 
loading series would, presumably, be of the same form as already shown 
for the subsonic case as given in eqs. (5).) The second difference 
between the sonic and subsonic cases lies in the ^orm of the kernel K, 
which is used for the calculation of the downwash factor F as defined 
in equation (27). The modification of the kernel K to the limiting case 
M = 1 has been performed in reference 2k, but the integration of the 
kernel with respect to the spanwise variable T) must still be performed 
and is presented in another part of this section. 

The present method is based on the usual assumptions of linearized 
theory and the use of this approximate theory may be open to question. 
However, for the very thin wings now being used on aircraft, it is 
felt that the first-order effects as given by linear theory constitute 
the major effects for unseparated flow and adequate solutions will be 
obtained. - 1 - 


1-Some conditions which must be met in the linearization of the 
governing equations in the neighborhood of M = 1 have been presented 
by Miles (ref. 36) and Mollo-Christensen (ref. 37)* As stated by Miles, 
it is required that the conditions 8, MB, KB, KM8 « 1 be satisfied 
and, in addition, that one or more of the following conditions be met: 

j M - 1 j » 5 2 / 5 k » b 2 ^ - » 8 1 / 5 

where 8 is the thickness ratio, M is the Mach number, k is the 
reduced frequency, and A is the aspect ratio. 
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For the case of a sonic trailing edge the form selected for the 
loading is 


f -I 


A 


—f - J ^QO + ^l 1 ! + a 02 T l 2 + * * *) + 


~ f 2e( a 10 + ‘^ll 1 ! + a 12 + • • •) + (i " f Z e )^ 2 ( a 20 + 


a^ 1 ! t ^pp^] + •••)+... 


•) 


(32) 


or 


£p _ 8 

q. t> 


where 

1 

f n(5- f Ze) = (6 - f j e ) 2 (n - 0, 1, 2, . . .) 

Anri 


n=0 m=0 


UJ 

/l. a rmi*> n "^2 e) ("H ) 


(33) 


%(n) = Ti m \/ s2 - ^ 


This series contains some aspects of the supersonic case as well 
as the subsonic. For instance, as In the supersonic case, the loading 
is not zero at the trailing edge; whereas, as in the subsonic case, the 
loading becomes infinite at the leading edge. 


Treatment of Integration in Chordwise Direction 

Dy following the scheme used for M < 1, a set of simultaneous 
equations for the chordwise replacement loads is obtained corresponding 
to equations (15) and (l6). 

In calculating the values of the replacement loads, it is convenient 
to choose the semichord at some spanwise station as a reference length, 
as was done for M < 1, and to perform a coordinate transformation such 
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that the origin lies on the wing midline. Ery such a transformation a 
new chordwise loading function f n (|) can he defined as 


f n (0 = (1 + 2 


(n — 0, 1, 2, . . .) 


where g. is now considered as a dimensionless variable based on a 
reference semichord and measured from the wing midline . In terms of 
this new variable g, the equations expressing the loading conditions 
(corresponding to eqs. ( 15 ) in the subsonic case) appear as 

j ol n-i 

£ G vn = / (1 + I) 2 d| (34) 

v=l J -1 

and for the downwash condition (corresponding to eqs. (l6) in the 
subsonic case) 


YZ G K = / (l + g) 2 K . dg (35) 

^ vngv J_ ± g| 

where the functions and designate 

v - M=1 

and 



and are the two-dimensional kernel as defined by equation (B23) of 
reference 24 and given in appendix B by equation (B9) . Unlike the 
two-dimensional kernel for the subsonic case, these functions are not 
tabulated for M = 1 and hence must be calculated. 

In expanded form and for the nth chordwise-loading term, the set of 
simultaneous equations is 
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Gin + G^ + G^n + G^ 



G jAl + G Jl2 + VS3 + Vl4 “£ X « 

G iAl + a sAs + Vfe + ^>*21 « 

Gin%L + g ^%2 + ^3x^33 + %n%4- = j/ 1 f n(S )%{ d £ 


A set of simultaneous equations as given by equation (36) must be formed 
for each value of n and solved for the. associated values of G. 


As in the case of subsonic flow, the integrands of the integrals of 
equation (36) become infinite at £ = x g- ® ie evaluation of these 

integrals is discussed in the latter part of appendix B. 

At sonic conditions, a disturbance cannot be propagated upstream, 
so that the labor of computing the loading functions is reduced as 
compared with the case for M < 1 since the downwash factors do not have 
to be determined and summed for conditions where the control point is 
upstream of the Integration area. 

As an example , if the loading functions are assumed to act at the 
l/8-, 3/8-, 5/8-', and 7/8-chord positions and the control points are 
located at the l/h-, l/2-, and 3 A -chord positions, then 


K 12 - K 13 “ K l4 “ *23 ~ ~ *34 " 0 

since the terms represent downwash factors for control points ahead 
of the loading stations. Thus, the set of equations becomes triangular 
and can be solved by successive substitutions. 


Treatment in Spanwise Direction 

The integration in the span direction is carried out in the same 
manner as was done for M < 1, that is, the wing is divided into many 
small segments, and, with the load Ap assumed to be constant across 



the segment, the integration is performed. The kernel for this case is given by equations (47) in 
reference 24. However, it is not possible to integrate this expression for the kernel in closed form 
and, as was done for the case -when M < 1, the kernel was expanded in powers of the frequency param- 
eter k and integrated term by term. The concept of the principal part of a finite integral was 
used. The form of the integrated series is given in the following equation: 




ik*o r 

1 1 


..3 r 


jjy, + . < y . dt] * jsL— ^ + - <y - 1) 3 ]} ♦ ^ (L o {rh - tH 

~ 2(1 - 7 - ~ + log - %-f^- - 7^1) - (y + Di°e(y + i) 2 + (y - i)iog(y - i) a + 
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■where the parameters Xq, y, k, and y are Identical to those 
employed in equation (31) . 


REMARKS PERTINENT TO THE SURFACE -LOADING METHOD 


In previous sections of this paper, some of the approximations 
involved in the present lifting-surface method have teen discussed. It 
has teen pointed out that certain arbitrary features arise in performing 
the numerical evaluations of the integrals, particularly -with regard to 
the chordvd.se integrations. The assumptions that have teen made in this 
connection include the replacement of the continuous chordwise loading 
by a number of Individual replacement loads, the use of the two- 
dimensional kernel function in evaluating these loads, and the choice 
of a set of load stations at which the replacement loads are assumed to 
act. In addition, a set of control points must be selected at which 
the governing downwash conditions are to be satisfied. In the actual 
application of the procedure these features give rise to certain problems 
which are discussed in the following paragraphs. 


Control-Point and Load-Station Location 

No attempt has been made in this paper to determine the most favor- 
able location of the control points or the optimum location for the sta- 
tions at which the replacement loads are assumed to act. For such a 
study, systematic tables of the downwash factor are desirable and when 
these tables become available, it will be much easier to evaluate this 
aspect of the problem. In the present study, such tables were not avail- 
able, and the positions of the control points and the means of distribu- 
ting the loading selected were the same as have been used for the steady 
case. It is felt that, at least for the simple rigid modes considered 
in the present study, the location of the control points is not critical. 
For more complicated modes of deformation other locations nnr) additional 
control points might be necessary. This question is worthy of further 
study. 


Effect of Taper 

In the calculation of the distribution of the loading across the 
chord, it is necessary to select a value of the reduced frequency k = — 
to be used in calculating the two-dimensional downwash factors pnfl 

K g £ in the Integrals in equations (l6). For a tapered wing, the problem 

arises as to what span position to use to obtain a reference chord for 
calculating k. 



36 


MCA TN 3694 


In relation to this effect of taper, as a part of the present 
investigation, two separate calculations were made for the case of a 
45° delta wing. In one case, the root chord was used as the reference 
chord, whereas, in the other, the chord at the midsemispan was used. 

It was found that the final results obtained for the pressure distribu- 
tion and, consequently, for the lift and moment from the two calculations 
were almost identical. Cm the basis of this one test case, it may be 
inferred that the location of the reference chord is not an .important 
factor. In the remainder of the delta-wing calculations presented in 
this report, the chord at midsemispan was used. 


DISCUSSION OF SOME APPLICATIONS OP THE METHOD 


This section is concerned with a discussion of the results of cal- 
culations which were based on the lifting-surface method discussed. 
Calculations have been made for both rectangular and triangular wings 
and comparisons of the results have been made with existing theoretical, 
and experimental results where possible. 


Rectangular Wing 

In order to furnish a basis for comparison of results of the method 
with existing theory, calculations have been performed at M = 0 for a 
rectangular wing with an aspect ratio of 2 pitching about the midchord 
for various values of the reduced frequency k. The results are compared 
with the results given by Lawrence and Gerber (ref. 29) and are shown 
plotted in figure 2. In figure 2(a) the magnitudes are plotted, and in 
figure 2(b) the corresponding phase angles are plotted. Excellent agree- 
ment is obtained for both the magnitude and the phase angle. Although 
both the present method and the method of Lawrence and Gerber are approx- 
imate, the good agreement between the results promotes a feeling of 
confidence in both methods. It should be noted that the two methods 
are not similar and contain entirely different approximations. 

To obtain some effects of Mach number, calculations were made for 
the same rectangular wing with aspect ratio of 2 pitching about the* 
midchord for a range’ of Mach numbers at a constant value of k = 0.22. 
Results are shown in figure 3- Included in the figure are the results 
of two-dimensional calculations . The magnitudes of the lift and moment 
are given in figure 3(a) and the corresponding phase angles are given 
in figure 3(b). Calculations for three-dimensional flow up to and 
including M = 1 were made by the use of the surface -loading method. 

The results at supersonic speeds were obtained from reference 12. The 
variation with Mach number for both the lift and moment is approximately 
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that which would he predicted by use of the factor & - up to 

M = O. 7 . Both the moment and lift increase at M = 1 and then drop 
off again at supersonic speed. Note that the two-dimensional lift 
coefficient has the same shape for the range up to M = O .7 and that 

A 1 

application of the aspect-ratio correction factor — — — = — would 

apply fairly well. No such simple factor exists for the phase angles, 
and it would not he possible to correct two-dimensional results for 
finite span effects. 


Triangular Wing 

Figures 4 and 5 present results for a delta wing with an aspect 
ratio of 4 oscillating in pitch about the midchord. In figure 4 is 
shown the lift and associated phase angles plotted against the reduced 
frequency k for M = 0. Results of the present analysis as shown by 
the solid line and those of Lawrence and Gerber (ref. 29) as shown by 
the dashed line are compared with some experimental results (indicated 
by the circles) obtained by Sumner A. Leadbetter and Sherman A. Clevenson 
at the Langley Aeronautical Laboratory. It is noteworthy that the results 
of the two methods agree rather well, even with respect to phase angle. 

In figure 5 a corresponding comparison for the moment and its asso- 
ciated phase angle (for the same wing and for the same conditions as in 
figure 4) is shown plotted against the reduced frequency k. Results of 
both analyses are in substantial agreement with respect to the magnitude 
of the moment. With regard to the phase angles, a significant difference 
between the results of the two theories occurs, although both theories 
indicate the same trend. The source of the difference can not be 
explained at the present time and will have to be resolved by further 
calculations and experiments. 


Results of Multiple -Line Method 

Figure 6 shows results of calculations for a rectangular wing with 
an aspect ratio of 2 at M = 0 based on the multiple-line method 
described in appendix. D. Results of the line method approach those of 
the surface-loading method when a fairly large number of lifting lines 
and control points are used. These results are discussed more fully 
in appendix D. 
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CORCLUDIHG REMARKS 


The purpose of this paper has been bo present and describe In 
detail a method for calculating the loading on a -wing of any practical 
plan form which is oscillating in a subsonic or sonic stream. The 
method is presented in general form and some results of application are 
discussed. A sample calculation for a specific case is presented in an 
appendix. The method may be used for calculating the loading on elastic 
wings as well as on rigid wings, ibis feature makes the method adaptable 
to flutter calcinations since it is possible to calculate the loading for 
the various modes us ually assumed for a normal type of modal flutter 
analysis. The method can also be applied in principle, at least, to the 
combined aerodynamic and structural problem in which the flutter char- 
acteristics are obtained directly by the use of structural and aerody- 
namic influence coefficients. 

The procedure used is based on linearized theory in which the usual 
assumptions of linearized flow, such as small thickness ratio of the 
wing, inviscid fluid, and so forth,- are necessary. 

The method has been found, to give good agreement with existing 
theory and experiment for low subsonic Mach numbers for both rectangular 
and triangular wings. Results for high subsonic Mach numbers and a Mach 
number of 1 fit in well with theoretical results for supersonic flow but 
require further verification by comparison with experiment. 

In addition to the presentation of the lifting-surface method, a 
multiple-line approach is included in an appendix. Results of the line 
method approach those of the surface -loading method when a fairly large 
number of lifting lines and control points are used. 

It is realized that several variants of the procedure may be made 
and may be desirable for routine or systematic calculations, particularly 
in view of the constantly increasing capabilities of automatic computing 
equipment. 


Langley Aeronautical Laboratory, 

Rational Advisory Committee for Aeronautics, 
Langley Field, Va., March 14, 1958. 
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APPENDIX A 


NUMERICAL EXAMPLE OF THE LIVING-SURFACE METHOD 
FOR A RECTANGULAR WING 
(A = 2; k = 0.22; M = 0.5) 


As an illustration of the calculation procedure of the lifting- 
surface method, the details of the procedure are presented in this 
appendix for a rectangular wing with A = 2 and k = 0.22, which is 
oscillating as a rigid wing about the midchord line in a stream flow 
of M = 0.5- For the case of subsonic flow the main equation to be 
dealt with is equation (30a) or (30b). The steps in performing the 
calculations are as follows: 

(l) The first step in the process is to divide the wing into a 
number of areas as shown in figure 1, to select the number and location 
of the control points, and, consequently, to determine the number of 
terms of the loading series (eqs. ( 5 )) which are to be retained. For 
the present case, the wing was divided into four equal chordwise areas 
for the calculation of the distribution of the load in the chord direc- 
tion. For the spanwise integration, the wing was divided into 21 seg- 
ments of which 19 are segments of span 2e and two (one at each tip) 
are of span e. The semispan e of a segment is thus equal to one- 
twentieth the wing semispan, (in all of the calculations this quan- 
tity e is used as a nondimens ionalizing factor so that the full wing 
span becomes 40 units. For a wing with A = 2, as considered in this 
example, the chord is therefore 20 units and extends from -10 at the 
leading edge to 10 at the trailing edge. The reduced -frequency param- 
eter k (based on the half -chord) in terms of e becomes k = ^ ^ 

b V 

and for this example k = 0.22 and k = 0.022.) 

Nine control points located at the l/k - , l/2-, and 3A~chord posi- 
tions and at three span stations tj = 0.2, 0.5, and 0.8 of the wing 
semispan were selected. The selection of nine control points determines 
the minimum number of terms of the loading series. There can be more 
control points than terms of the loading series; however, a method of 
least squares would have to be employed for the final solution. 

In the present example, nine terms were retained in the loading 
series. Three chordwise terms, each modified by three spanwise terms, 
were used so that the series contains nine unknown coefficients a Tim and 

appears as 



_ 8 / 2 
<1 b V 


8- - TJ 


cot — fa + a + a , + sin 9 (a + a n‘ 

2 ^00 02 Ok J ^10 12 


•1 


+ a . tH 4- din 20fa + a t) 2 + a tA, 
Ik [ 20 22 2k 


(Al) 


(A2) 


It is only necessary to consider the even power terms in q since the loading Is assumed to be sym- 
metrical in the span direction about the mldspsn position. 

(2) The next step consists of calculating the replacement loads G from equations (17), (l8), 
a n d (19) . For the present case, there are three (n = 0, 1, or 2) sets of four simultaneous equations 
as follows: 

°ln + a 2n + + °4n “ J^ f n^ ^ 

(-2.83500 - 0.450001)0^ + (2.25500 - 0.476801)0^ + (0.60680 - 0.296701)0^ + (0.29180 - 0.220661)0^ - f * n (i)S^ 46 

(-1.17430 - 0.217601)0^ + (-2.83500 - 0.450001)0^+ (2.23300 - 0.476801)0^ + (0.60680 - 0.296701)0^ - J f^l)^ di 

(-0.85000 - 0.090371)0^ + (-1.17430 - 0.217601)0^ + (-2.83500 - 0.45000i)G 3n + (2. 25500 - 0.47680l)a^ n - J f n (s)% { ^ 

The coefficients of the replacement loads G are values of the two-dimensional kernel and 

have been obtained directly from the table of reference 31. The kernel ifon Is a function of M, and 
Z » k(x - |). Consider as an example the center of area at l/8 chord ■= anji a corr ^ ro1 point at 

l/k chord (x = -5), For k » 0.022, Z » (0.022) ^5 + ~) = O.O55 and for M = 0.5, the downwash 
factor K u - -2 .83500 - 0.450001 may be read from the table of reference 34. Similarly, for the 
same values of E and M and for the center of area at 7/8 chord «= and the control point at 

midchord (x *= 0), Z = 0.022 |o - = -O.165 and, therefore, K^ 2 “ 0.60680 - 296701. 


-r 

o 


The integrals on the right of equation (A2) axe given aB follows for the three sets of equations 
correspo ndin g to n = 0, 1, and 2. (See appendix B for discussion of method of integration.) 
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For 


For 


n = 0, f 0 (|) = cot — and 



n = 


cot — sin 0 d0 = jt 
2 

cot | K-^ sin 0 d0 = - 2.88354 - 1-341251 

cot | K 2| sin 0 d0 = -2-97504 - O.96137I 

0 

cot g sin 0 d0 = -3.03074 - 0-573201 

1, = sin 0 and 


j 



0.58912 - 0.63606i 

-0.45095 - 0.670321 
-1.48547 - 0.553661 




> 


4l 


(A3) 


(A4) 



42 


NACA TN 3694- 


For n = 2, f 2 (|) = sin 26 and 



sin 26 sin 8 d8 = 0 


sin 28 sin 8 





= •-1.03882 - 0.20631i 


> 



sin 28 sin 8 d0 = -2.02994 + 0.04l45i 
sin 28 sin 8 K_ d8 = -1.02892 + 0. 288611 


(A5) 


The Values of G as computed from the three sets of simultaneous 
equations formed "by substituting equations (A3), (A4), and (A5) into 
equation (A2) are 


= 1.72184 + 0.00224i" 

G 20 = °* 73737 ' 0*011301 
Gjq = 0.44052 - 0. 00484 i 
% 0 = 0.24186 -f 0.0335891 

G^ = 0.30752 + 0.022411 
G 21 58 0.47670 + 0.005821 
G^ = 0.47702 - 0.005831 ' 

= 0.30956 - 0.022401 

g 12 = 0*^5744 + 0.011031 
g 22 = 0.23804 - 0.010551 
G32 * -^O. 23 753 - 0.010731 
G 1^2 = -0*45795 + 0. 010241^ 

When the values of G are obtained, the loading function AP v (p) 

appearing in equation (30a) and defined by equation (26) may be written 
as 
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AP v (n) 



+ a o # 2 





+ G 


v2 



+ a 2i . 2 



(AT) 


where v =1, 2, 3, or ^ corresponding to the k areas in the chord 
direction . 

Examination of this equation shows that certain products of 

G vn^ n/i " 7 p 2 are needed. These products correspond to the products 

®m.(^) G vn sduation (30h). It has been found convenient to arrange 

these products in a certain form for later calculations . This form is 
shown in table H where values of the products are given for the complete 
system of replacement loads. 

(3) The downwash factors F must now be determ in ed. These factors 
may be calculated by the use of equation (31) , which was used for this 
example. With the adaptation of high-speed computing machines to the 
problem, the kernel in the form of equation (20) of reference 2^ may be 
numerically integrated without, perhaps, the frequency and Mach number 
restriction of the series. However done, the distance between the centers 
of areas and the control points, x Q and y', must first be determined. 
For the present case, the y’ distances for control point 1 and the 
first chordwise area corresponding to g = 1, for example, are 


u 

0 

0.1 

0.2 

0.3 

0.1)- 

0.5 

0.6 

0.7 

0.8 

0.9 


D 

2 

0 

2 

k 

6 

8 

n 

12 

14 

m 

■ 

6 

8 

10 

12 

lA 

16 

H 

20 

22 


The distance x = 2.5 is constant since the wing is rectangular. 

Since the loading in the example is symmetrical about the midspan, 
the loading factors which multiply F in equation (30b) are the same 
for equidistance on each side of the midspan. Therefore, the values 
of F which are located at equal distances from the mid span may be 
computed and added before multiplication by the loading; consequently, 
distances y' R and y' L have been given In the example where y' R 

refers to the distance from a control point to the centers of the areas 
to the right of the wing midspan and y'-^ refers to distances to areas 

to the left of the midspan. The final form of the F factors is given 
in table TTT where, for each entry, the two values of F for yi^ and 

y’k have been added. 




















I 


( 4 ) The next step la the multiplication of the loading factors of table XI by the downwash 

factors F of table HI in accordance with equation ( 50 b) which results in a 9 "by 9 matrix of complex 
elements. This is given by the left side of the following equation and the downwash conditions to be 

discussed subsequently are shown on the right of the equation: 

-\ 

(l.omj + 0.1&L56i)aQQ + (-0.05279 + 0 , 09569 l)a 10 + ( 0 . 35112 + O.Cfcn^Jago + (- 0.14305 + O.o43071)ttt}2 + (- 0*05511 + 0.014651)832 + 

(- 0.02878 + 0 . 007161 )^ + (- 0 . 0 Tl (*9 + 0 . 018031 )b^ + (- 0.02325 + 0 . 005881 ) 0 ^ + (- 0.01569 + 0.002621)0^ - -( 0.025 - 0 . 002751 ) 2 * 

(1.20324 + 0 . 02339 i)a 00 + ( 0.27977 + 0.086131)1^ + (0.60906 - 0.031*21)820'+ (-0. 2199*1 + 0.061731)002 + (-0.06529 + 0.022081)833 + 

(-0,03123 + 0.009211)822 + (-O.IO617 + 0.027581)8^ + (-0.03939 + 0.0092$H)*j^ + (-0.0190+ + O.OOt^i)^ a -(0.025 + 01)2* 

( 1.28988 - 0.123231)8^, + ( 0.60852 + 0 . 031521 ) 8 ^, + ( 0.32497 - 0 . 104811 ) 8 ^ + (-O.26974 + 0.0^7921)8^ + (- 0.11441 + O.OSSAi)^ + 

(-0.02668 + 0.011771)822+ (-0. 12604 + 0.04013l)a<*+ (-0.05509 + 0.017*71)8^ + (- 0 . 00.473 + 0.006151)82^ - -(0,025 + 0.002751)2* 

(0.90622 + 0. 110761) + (0.01496 + 0.066721)8^ + (0.27334 + 0 . 022491 ) 8 ^, + (0.99139 + 0 . 001481 ) 8 ^ + (0.11217 + 0.022721)8^ + 

( 0.25999 + 0.003931)822 + (0.69160 - 0 . 007151 ) 8 ^ + ( 0.08428 + 0.002261)8^ + (0.17638 + 0.00n6±)e^ - -(0.025 - 0.00275l)2rt 

(1.31514 - 0.019621)800 + ( 0.29357 + 0.052541)830 + (0.47578 - 0 . 034101)820 + ( 1.24921 - 0.1301015802 + ( 0.37621 - O.OOlSEi)*^ + l ^ 

{ 0.41833 - 0.042711)822 + (0.86358 - 0,096661)*^ + (0.26493 - 0 . 004901 ) 8 ^ + (0.28075 - 0.029771)8^ - -( 0.025 + 01 ) 2 * 

(i.© 34 l - 0 . 156241 ) 8 ^, + (0.56690 + 0.000261)8^ + ( 0.26664 - 0.068061)0^ + (1A2203 - o.2ffl.d9l)8 02 + (0.63972 - 0.060631)832 + 

(0.25104 - 0.068381)822 + (0.98606 - 0.201641)8^ + ( 0.44244 - 0.045991)8^ + (0.17013 - 0 . 06026 l)a^ - -(0.025 + 0.002751)2* 

{1.01349 + 0.157721)8^ + (-0.05208 + 0.065731)83^, + ( 0.30911 + 0.035201)*2q + (O.25606 + 0 . 059421 )» TO + (- 0.00453 + 0.019791)832 + 

(0.07696 + 0 . 006151 ) 8 ^ + (- 0.04884 + 0 . 023251 ) 8 ^ + (- 0.02059 + 0 . 007971 ) 8 ^ + (- 0 . 0069 + + O.OO 440 i)ag, ° -( 0.025 - 0 . 002751 ) 2 * 

(1.15764 + 0.009751)8^ + (0.27905 + 0 . 0748 u)e 10 + (0.56079 - 0.0538U)b 20 + (0 . 29204 - O.OOlSaOa^ + ( 0.07234 + 0 . 036611 ) 8 ^ + 

(0.13788 - 0.006791)822 + (-0.09262 + 0.028611)8^ + (-0.03666 + 0.011331)8^ + (-O.OC479 + .0.003931)8^ - -(0.029 + 01)2* 

(1.24922 - 0.151261)8^, + (0.58649 + o.o2i&7i)a 10 + (o.soero - 0.090261)820 + (0.31639 - 0.037201)002 + ( 0.14026 + 0.003211)832 + 

( 0.07499 - 0.024521)822 + (- 0 . 1159 + + 0 . 040071 ) 8 ^ + (- 0.04055 + 0 . 016031 ) 8 ^ + (-O.OC 0 OO + 0 . 004021 ) 8 ^ - -(0.025 + 0 . 002751 ) 2 * 

( 5 ) She bou ndar y conditions w at the various control points are determined from the motion of 
the wing. For the present case, oscillations about the midchord line of a nondeforming wing ere 
considered. 
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Since f(x,y) 


a(x - a) 


and, according to equation (3a), 


v(x,y) 


v 5f (x,y) + ^ ( x >y) 

bx 8t 


then 


v(x,y) = - [Va + a(x - a)] 
and for harmonic motion 


w(x,y) 

Va 


= - [l + ik(x 



(A9) 


It is convenient to divide the simultaneous equations by the factor 
4o/2rf which appears in equation (30b). Therefore, the downwash condi- 
tions for rotation about midchord (a = 0) are 


w ( x l’ y l) 


w(x7,y T ) 

Va 

Va 

Va 

w(x 2 ,y 2 ) 

u(x 5 ,y 5 ) 

v (*a’*a) 

Va 

Va 

Va 

v ( x 3’ y 3) . 

_ w (V y 6) 

w ( x 9^9) 

Va 

Va 

Va 


-(0.025 - 0.00275i)2rt 


-(0.025 + 0i)2it 


f (A10) 


-(0.025 + 0.00275i)2rt 


where x^ and y^ are the coordinates of the control points. 

(6) By the use of the boundary condition w/Va, the set of simul- 
taneous equations (eqs. (A8)) may be solved for a^. For the example 

considered the values of the a ^ m are 


Sqq = - 0.150858 + 0.033^381" 
a 1Q = 0. 0i)-5628 - 0.1291621 
a 2Q = -0.000862 - 0.0083371 

a 02 = -0.034965 + 0.0071071 

a l2 = °-°5T334 - 0.0316151 r 
&22 = 0.014221 - 0.004 922i 

a o4 = -0.025231 + 0.0038601 

a^ = 0.058109 - 0.006l98i 
= 0.045049 - 0.007729i 


(All) 
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( 7 ) These coefficients may be used in equations ( 5 ) to determine 
the pressure distribution and. In turn, to obtain section forces, total 
forces, and so forth. The total lift and moment, for example, are 
given by equations (8) and ( 9 ) and for the sample case are 


L,aj = 2.632 


% 


cq 


= 1 . 59 ^ 



19^ A3° 
3 ^ 9 - 14 ° 


(A12) 
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APPENDIX B 

TREATMENT OF CERTAIN INTEGRALS WHICH CONTAIN SINGULARITIES 
THAT ARISE IN THE CHDEDWISE INTEGRATION 


In "the treatment of the chord-wise Integration in the surface- 
loading method, certain integrals arise -which contain singularities 
(eqs. (l6)) and which must he given special treatment. In the deter- 
mination of the coefficients appearing in these equations, it is 

necessary to evaluate integrals of the form 

^ = jL* as (m) 


where f n (£) is the nth typical chordwise pressure mode in the 

series_expression for the continuous pressure distribution and 
Kg^ = KjM,k(xg-£^] is the two-dimensional kernel function. In this 

equation, the kernel function becomes infinite at £ = x g so that the 

integrand is singular. In the following sections, the methods of eval- 
uating the integrals and taking care of the singularities are discussed. 
First the case of ■ subsonic flow and then the case of sonic flow is 
discussed. 


Case of Subsonic Flow 

For the case of subsonic flow, the integrals occurring in equa- 
tions (l6a) and (l6b) are 



cot 7 T sin 0 K t d9 
2 gi 


(B2a) 


I = 
n 


'0 


sin n0 sin 0 K . d0 
Si 


(n = 1, 2, . . .)(B2b) 


at 


Since the kernel in the integrands of Iq and 1^ becomes singular 
£ = Xg, it is necessary to separate the kernel into singular and 
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nonsingular parts. The singular part can be integrated in closed 
analytic form. The nonsingular part can be handled routinely and 
accurately by numerical means. The separation of the singularity has 
been accomplished by Schwarz (ref. 3*0 in the following manner: 


+ K^Z) (B3) 


K(M,Z) = ^ F(M) + iG(M)log Zl 
Z G 


where 


Z = k(x g -|) 

The singularities now appear only in the quantity in brackets, and the 
nonsingular part K-. has been tabulated by Schwarz (ref. 3*0 • Substi- 
tuting this expression into the integral of equation (B2a) gives 
( after setting | = -cos 0 and x g = -cos '01 


/ 

Jo 


cot — sin 0 K(M 
2 


,Z)d0 = r* + cos 0)K L (M,Z)d0 + 

+ COS + iG(M)loge|zj] d0 (b4-) 


The first integral on the right-hand side of equation (B*0 must be 
evaluated numerically. The second integral may be found analytically 
and has the following value: 


/> 

J 0 


+ cos 0) 


F(M) 


iG(M)logg |z| d0 = ^ F(M) + iG(M) loge 2 - 


rt cos 9 + it 


l°ge k] 


(B5) 


The treatment of the integral of equation (B2b) for n = 1 proceeds 
in a similar manner to yield 
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pit pit 

/ sin 2 0 K(M,Z)d0 = / sin 2 © K-j(M,Z)a0 - ^ f(M) cos 0 o + 

J 0 JO k 

G(M) ^-log e 2 + i- cos 20 o + log e kj (B6) 


Correspondingly, for n = 2, 


prt pit 

J q s in 20 sin 0 K(M,Z)d0 = J sin 20 sin 0 K x (M,Z)d0 - £ F(M) cos 20 o 

^ G(M) ^ - y - 9 - ° - cos 0 o j (B7) 


Case of Sonic Flow 


For the case of sonic flow, the form, of the series expression for 
the chordwise loading and also the form of the kernel function differ 
from those of the subsonic case. 

For M = 1, the set of equations to be solved for the chordwise 
loads Gyn is given by equations (34) and (35)* In equation (35 ) , 
integrals of the following form appear: 


^ = 



n_ 2 

+ 0 ^gi 


ag 


(n = 0, 1, 2 . . .) (B8) 


At M ='l, the two-dimensional kernel function appearing in these 
equations is given by equation (B23) of reference 24 and may be written 
as 



where 
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and where 




are Fresnel integrals defined by- 


cos t%t 
2 


and 


S(a) 




dt 


As in the case of subsonic flow, the kernel function can be sepa- 
rated into singular and nonsingular parts. For this purpose, equa- 
tion (B9) may be rewritten as 



The integral 


In ( e d* (B8)) can then be rewritten as 



(1 + I) 





(B13) 
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[Che first integral on the right of equation (B13) is nonsingular 
and can he evaluated numerically. The second integral contains a 
singularity and can he integrated analytically. 

For n = 0, the second integral on the right of equation (B13) 'is 
found to he 


f* g 


as 


/ - ----- = 

j_l \j 1 + S \/(l + S)(x g - I) 


(Bib) 


where 


c - fx/f (1 + 


(B15) 


Correspondingly, for n = 1, the integral becomes 


°£ g ^=3 at=0 k + 1 )f 


(Bl6) 


For n = 2, there is obtained 


£±±^l u = c{xs + rf 


n / 


Xg - S 


(B17) 
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APPENDIX C 

CALCULATION OF MOMENT ON DELTA WING 


As mentioned in the text, the expression for the moment depends on 
the particular plan form. As an indication of the general method, the 
specific case of moment on a 45° delta -wing is given in this appendix. 
The moment is calculated about a line | = a as shown in the following 
sketch: 





The expression for the moment as given by equation ( 7 b) may be 
written as 


M a = qSbC^ a =JJ (6 - a)Ap d| 


dn 


(ci) 


When equation (k) is used (for a ^5° delta wing) the distance of 
an arbitrary point (ti,£) from the line £ = a in terms of the vari- 
able 9 is 
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a = s 

§ tan A - 

f- 

- ^ tan A|cos 9 

*o~ 
- a — 


2 

V 8 

2 . / 

s_ 


(C2) 


Substituting Ap (from eq. (5a)) and equation (C2) into equa- 
tion (Cl) and noting that 


= s^a 1 tan A^sin 0 d0 = sb(T])sin 0 d0 


results in the following: 


Ma'SB^pv - I 2 =ot f(a 00 + s^ 2 + + 

ato 6 (a 1Q + a^n 2 + + sin 26 (a^ + a^r, 2 + a^) 


tan A - 


tan a) cos 0 - a — 
K* 2. J s 


sin 0 d0 dT] 


(C3) 


In this example, three chordwise and three spanwise terms have been 
retained in the series for the loading Ap. 

Performing the indicated integration in equation (C3) results in the 
following equation for the moment on a 45° delta wing ( by mpTH-ng use of 

the fact that tan A = 2 ^ 
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APPENDIX D 

DESCRIPTION OF A MULTIPLE-LINE METHOD 


Since a lifting-line method of handling the integral equation 
(eq.. (l)) may often be simpler than a lifting- surface procedure, it was 
considered interesting to investigate such a method and to compare 
results with those obtained by the surface approach. The purpose of 
this appendix is to present a lifting-line method based on a steady- 
state procedure developed by Schlichting and Ka hl ert (ref. 26) and to 
give some results obtained by it. 


Description of the Method 

In the Schlichting procedure for steady flow the lifting surface is 
divided into n spanwise strips of chord c‘ = c/n, where c is the 
total chord. A lifting line is placed at the l/4~ chord position (c 1 A) 
of each strip and the downwash is satisfied at control points on the 
3 A -chord position of each strip. This 1 A -chord and 3 J^- -chord. location, 
respectively, of the lines and control points is an essential element of 
the procedure . It is a known fact that for a two-dimensional flat plate 
in steady incompressible flow the exact value of the lift can be obtained 
by placing a single lifting line at the lA-chord position and satisfying 
the downwash at the 3 A -chord position. Schlichting and Kahlert made use 
of this fact in developing a procedure for finite swept wings . 

In adapting Schlichting ' s method to the oscillatory case the same 
placement of lifting lines and control points has been used and is shown 
in the following sketch, where, for example, the wing has been divided 
into two spanwise strips: 
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In this example two lifting lines and 1 2 are considered 

that lie, respectively, at l/8 and 5/8 of the total chord. Two control 
points are shown, by Xp at 3/8c and X£ at j/8c. In the sketch 4> is 

an angular spanwise coordinate related to a nond imens ional variable tj 
(referred to the semispan) by t) = -cos $ . 

In the original form of the integral, equation given by equation (l), 
the continuous pressure distribution Ap is replaced by the sum of 
loadings on the individual lines. It is assumed that the loading on the 
line Z n can be represented by a series of the form 

Apn = pV 2 s ^a^ sin <5 + a ^ sin 3® + a^ sin 5$ + . . A (Dl) 


which contains only one variable, the spanwise coordinate <p. S inc e a 
separate loading function of this form is written for each lifting line, 
the double integration of equation (l) is reduced to a sum of single 
integrals. The integral equation then appears as 


v(x,y) 

V 


8itqs^ 


t/. 


APn 


span 



(D2) 


where x v is the chordwise coordinate of the nth lifting line, and where 
K, the nondimens ional form of the kernel function defined by equation (29), 
has been employed . 

When Apn in equation (D2) is replaced by the series expression in 
equation (Dl), the integral equation is reduced to a summation of definite 
integrals multiplied by the coefficients a^, a Q ^, and so forth, as 

follows : 



where 


x 


X - Xy 


y Q = y - n 
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This equation sums the downwash at a particular control point due 
to all the lifting lines and the various types of loading, sin $, 
sin 3 ®, and so forth, on each line. The specific problem is the deter- 


mination of the coefficients a 


hi* n3 J 


and so forth. In order to 


obtain the coefficients a, equation (D3) must be -written for each of a 
number of control points. This leads to a set of simultaneous equations 
which can be solved for the values of a. At each control point the 


function 


v(x,y) 

V 


is determined from the motion of the wing in the same 


manner as in the surface-loading method. 


Once the values of the coefficients a^i , and so forth, have 

been found, they can be used in equation (Dl) to define the pressure 
distribution on any lifting line and in turn to give various force or 
moment coefficients. As an example, the total lift on a rectangular 
wing can be obtained from the relation 



(*) 


where 1^ is the lift on the nth lifting line. The total pitching 
moment about an axis x = a for a rectangular wing is given by 


C M = 


Ma_ 

qSb 



(D5) 


where x n is the chordwise coordinate of the nth lifting line. 

The' main computational problem in the procedure is the evaluation 
of the spanwise Integrals of equation (D3), and the handling of these 
integrals is discussed in the next section. 


Evaluation of the Spanwise Integrals 

The spanwise integrations in equation (D3) are performed numer- 
ically and, since the kernel function has not been tabulated, it is 
necessary to m ake use of the series form of the kernel given by 
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jd 


equations (31) and (5*0 of reference 2k. In. this form the nondimensional 
kernel function K appears as 



where all distances are referred to the semi span s. 

For the integrals of equation (D3) which relate to lifting lines 
behind the control point (x Q < , the products of the loading modes and 

the kernel K as defined by equation (D4) can be evaluated at a number 
of values of T) and the numerical integrations can be readily performed. 
Men the lifting line lies ahead of the control point (xq > 0) , however, 
the kernel becomes inf inte at f\ = y, and the singularities must be 
carefully treated. 

For the case x Q > 0 with q = y, singularities arise in the terms 



Extraction of the singularities .- The s ingular ities of the terms A 
and B can be extracted and, when combined with the loading terms sin $, 
sin yb, and so forth, are integrable in closed form. 
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For this purpose, term A can be rewritten as 



•where the quantity in brackets is finite for all values of r\ and where 



Tbe remaining term on the right of equation (d 8) contains the singularity 
and will be handled ana lytic ally . 


Ihe singularity of term B can also be isolated by writing 



(DIO) 


where the bracketed quantity is finite for all values of tj and where the 
remaining term can be handled analytically. 

Special form of the kernel for the' case > 0.- The expressions 

for the terms A and B obtained in equations (D8) and (DIO) can be 
used to define a special form of the kernel function for use in the case 
x o For this case the kernel can be written as 


K (M’ toI °’ k r O )x o>0 - K(* 0 >°,y 0 ) + K' (*o> 0 ,y o ) 


(Dll) 
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The term K^x o >0,y o ^ is nonsingul a r an d is defined by 


K(v°.y 0 ) - k(m.^ 0 ,^ 0 ) 


+ e'^o < 


k£ 

2 


■ 

io S - * 0 2 ) + 


(D 12 ) 


where K^M,k,Xo,y 0 j is the form of the kernel defined by equation (D8) . 

The second term on the right-hand side of equation (Dll) contains 
the singularities and is defined by 



-ikxp / 2 

\y 0 2 




(D13) 


The form of the kernel given by equation (Dll) is used in evaluating the 
integrals of equation (D3) for Xq > 0. 

Performance of the spanvrise integrations .- It is recalled that for 
the case x Q < 0, the integrals of equation (D3) can be readily evaluated 

numerically by making use of the form of the kernel given by equation (D8) . 
For the case x Q > 0, the integrals of equation (3) can be handled by- 
using the special form of the kernel defined by equation (Dll) . With this 
expression, a typical integral of equation (D3) for Xq > 0 can be written 

as 


sin m$ I^MjkxQ,}^^^ dq = J sin K^X^y^dq + 


fl 


sin m3) K' ^x 0 X),y 0 ^ dT) (D14-) 
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where the functions K^>0,y o ) ana K’ ^>0,y o J are defined, 

respectively, by equations (D12) and (D13) . Hie first integral on the 
right of equation (Dl^) can he readily evaluated numerically. Hie sec- 
ond integral contains the singularities and must be evaluated analyti- 
cally. Its values have been determined for two pressure modes, 

sin G> = sjl - T) 2 and sin 3$ = (^ 2 - l) \/l - tj 2 and are 


- -e' 1 ^ + 


J x/l - n 2 S’ (x o >0,y o )dTi 

k 2 \Jl - y 2 |(y + l)log(y + 1) - (y - l)log(y - 1) “ 2 J + 

tr (ii - n 2 - n / 1 - y^)i°g y G 2 dT i 


(D15) 


and 


sin 3 < I , 0 

— OTt — j 

sin «j> 0 


^ /l - t! 2 ^ 2 - l)K'(x 0 X),y 0 jdTi = -e- ikx o j 
k 2 ^1 - y^y 2 - lj £(y + l)log(y + 1) - (y - l)log(y - 1)- 2 J + 

Ji - n 2 (^ 2 - 1) -x/ 1 - y 2 ^ - x ) 


£ 

2 


log y D dii 


(Dl6) 


where 


* -1 
$ 0 = cos y 


For the special case of a control point at y = 0 and x Q > 0, equa- 
tion (D15) becomes 



K'(x o >0,-T])dTi = 


_ e -ikXo _ at _ + * log 2) 


(D17) 
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and equation (Dl6) becomes 


K - x ) 

For the case of y f 0 and 
of equations (D15) and (Dl6) 


K' (x 0 X»,-T]Jdii 


= -e" 1 ^ 


& + f (?) 


(D18) 


Xq > 0, the integrals remaining on the right 
are finite and can be evaluated numerically. 


Once the integrals of equation (D3) have been evaluated, the 
simultaneous equations in terms of -the unknown coefficients a m , a^, 

and so forth, can be formed and solved in the manner indicated earlier. 
In the next section application or the multiple-line method is made and 
results are compared with those of the surface -loading method. 


Application of Method and Discussion of Results 

The multiple-line method just discussed has been applied to the 
same rectangular wing with aspect ratio of 2 which was treated by the 
surface -loading method. Calculations have been made for several values 
of the reduced-frequency parameter k at M = 0 for the wing oscillating 
in pitch about its midchord. Two sets of calculations were made, one 
using two lifting lines and two control points and the second using four 
lifting lines and eight control points. In the two sets of calculations, 
the control points were located chordwise in accordance with the 3 A -chord 
concept discussed previously. With regard to spanwlse location, in the 
first set of calculations, the two control points were placed at the 
center of the wing; in the second set of calculations, four control points 
were at the center of the wing and four at 0.866 Bemispan. Results of the 
calculations are shown in figure 6 as the lift and moment coefficients and 
their associated phase angles plotted against k. The results of the 
surface -loading method are included for conqoarison. 

For the calculations of the lift based on two lines and two control 
points, the lift magnitude agrees well with the surface -loading results 
only at the lower frequencies. For the four-line, eight-point solution, 
results for the magnitude are approaching the results of the surface- 
loading method. Lift phase angles are in fairly good agreement for all 
sets of calculations ♦ 

With regard to the magnitude and phase angle of the moment, results 
of the line approach with four lines and eight control points are in fairly 
good agreement with those of the surf ace-loading method. There is signif- 
icant improvement over the results of the two-line, two-control point 
calculations . 

In general, it appears that in order to obtain accurate results with 
the multiple-line approach, a fairly large number of lifting lines and 
control points must be used. 
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TABLE I.- VALUES OF 
[n =* 0, 1, 2; v - 1, 2, 3, 5 

(a) M ■** 0 


B 

So 

G 20 

°30 

° 4 o 

HOOOOOOOOOOOO 

8 VS 8 ?^&>S^ 8 gaS 

1.718194 + 0.0006161 
1.719093 + 0.0014791 
1.720391 + 0.0018451 
1.724451 + 0.0006701 
1.728509 - 0.0031311 
1.731756 - 0. 0090461 
1.733788 - 0.0166151 
1-734444 - 0.0254521 
1,733669 - 0.0352721 
1.731489 - 0.0458541 
1.727947 - 0 . 05705 U 
1.723123 - 0.0687551 

0.736332 - O.OOO638I 
0.736393 - 0.0019901 
0.736438 - 0.0034461 
0.736223 - 0.0075211 

0.735275 - 0.0121351 
0.733421 - 0.0171551 
0.730576 - 0.0224531 
0.726706 - 0.0279251 
0.721815 - 0.0334841 
0.715925 - 0.0390491 
0.709075 - 0.0445621 

0.701300 - 0.0499821 

0.441765 - 0 . 0004121 
0.441523 - 0.0011621 
0.441159 - 0.0018161 
0.439865 - 0 . 00304 ® 
0.438213 - 0.0®732i 
0.436363 - 0.0039831 
0.434382 - 0.0038641 
0.432304 - 0.0034171 
0.430143 - 0.0026931 
0.427886 - 0.0017121 
0.425536 - 0.0004901 
0.423081 + 0.0009681 

0.245302 + 0.0004341 
0.244586 + 0.0016731 
0 . 2436 ® + 0 . 0034 m 
0.241054 + 0.0098911 
0.839596 + 0.0189981 
0. 240053 + 0.0301841 
0.242847 + 0.0429321 
0 . 248:39 + 0.0567941 
0.255966 + 0.0714491 
0.2662® + 0.0866151 
0.279033 + 0.1021031 
0.294089 + 0.1177691 

k 

^Ll 

G 21 

Si 

Si 

0 

0.02 

0.06 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.80 

0.90 

1.00 

0.306807 + 0.0015701 
0.306843 + 0.0046951 
0.306916 + 0.0077821 
0.307109 + 0.0152721 
0.307572 + 0.0224051 
0.307077 + 0.0291521 
0.306637 + 0.0355661 
0.305858 + 0.0416671 
0.304718 + 0.0474741 
0.303222 + 0.0530111 
0.301382 + 0.0582871 
0.299211 + 0.0633091 

0.478591 + o.ooo 4 ®i 
0.478544 + 0.0012041 

0.478430 + 0.0019981 
0.477887 + 0.0039321 
0.476691 + 0.0058071 
0 . 4755 ® + 0.0075621 
0.47363 9 + 0.0092631 
0.471315 + 0.0109021 
0.468535 + 0.0124811 
0.465307 + 0 . 0 l 4 oo 8 i 
0.461644 + 0.0154881 
0.457560 + 0.0169221 

0.478591 - 0.0004031 
0.478548 - 0.0012041 
0.478448 - 0.0019981 
0.478021 - 0.0039281 
0.477265 - 0.0057671 
0.476441 , 0.0074961 
0.475337 - 0. 0091141 
0.474033 - 0.0106131 
0.472541 - 0.0119851 
0.470852 - 0.0132241 
0.468970 - 0.0143251 
0.466894 - 0.0152781 

0.306807 - 0.0015701 
0.306861 - '0.0046951 
0.307002 - 0.0077821 
0.307779 - 0.0152761 
0.309268 - 0.0224451 
0.311775 - 0.02921B1 
0.315183 - 0.0357151 
0.319590 - 0.0419561 
0.325002 - 0.0479701 
0.331415 - 0.®37951 
0.338800 - 0 .® 9450 i 
0.347131 - 0.0649551 

k 

°12 

G 22 

G32 

°42 

HOOOOOOOOOOOO 

8883 g"j$ SS 8 SR 8 

0.460172 + 0.0014461 
0.460089 + 0.0021791 
0.459935 + 0.0036651 
0.459303 + 0.0074771 
0.458392 + 0.0113941 

0.457251 + 0.0153721 
0.455894 + 0.0193751 
0.454331 + 0.0233721 
0.452569 + 0.0273351 
0.450610 + 0.0312441 
0. 448453 + 0.0350831 
0.446108 + 0.0388351 

0.239286 - 0.0007151 
0.239261 - 0.0021401 
0.239193 - 0.0035591 
0.238877 - 0.0070741 
0.238383 - 0.0105401 
0.237728 - 0.0139541 
0.236911 - 0.0173091 
0.255937 - 0. 0206011 
0.234811 - 0.0238271 
0.233530 - 0.0269761 
0.232102 - 0.0300471 
0.230529 - 0.033® 01 

- 0 . 2392 86 - 0.0007171 

- 0.239261 - 0.0021541 
-O.239193 - 0.0035951 
- 0.238872 - 0.0072081 
-O.238368 - 0.0108241 
-O.237693 - 0.0144221 
- 0.236849 - 0.0179921 
-0.235838 - 0.0215181 
- 0.234665 - 0.0249811 
- 0.233331 - 0.0283781 
-0.231838 - 0.® 16961 
- 0.230195 - 0.0349231 

-0.460172 + 0.0007121 
-0.460089 + 0.0021151 
- 0.459955 + 0.0034891 
-0.459308 + o.oo 68 ®i 
- 0.458407 + 0.0099701 
- 0.457286 + 0. 0130041 
- 0.455956 + 0.0159261 
- 0.454430 + 0.0187471 
-O.452715 + 0.0214731 
-0.450809 + 0.0241101 
-0.448717 + 0.0266601 
- 0.446442 + 0.0290181 
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(TABLE I.- VAHJES OP 0 - Continued 

vn 

(b) K «■ 0.5 


\ °20 °30 


0000000000000 

1.718217 + 0.0008341 
1.718621 + 0.0015321 
1.719202 + 0.0020881 
1.719851 + 0.0025051 
1.720582 + 0.0027961 
1.721371 + 0.0023411 
1.722223 + 0.0029361 
1.723147 + 0.0027951 
1.723962 + 0.0025641 
1.724866 +.0.0021771 
1.728804 l 0.0013381 
1.731636 - 0.0071791 
1.732840 - 0.0146691 

OOOOOOOOOOOOO 

1 1 1 1 1 1 1 1 1 1 1 1 1 
OOOOOOOOOOOOO 

lifistitsspt 

0.441762 - 0.0005511 
0.441645 - 0.0010751 
0.44 i 498 - 0.00157H 
0.441332 - 0.0020371 
0.441132 - 0.00*791 
0.440915 - 0.0028961 
0.440675 - 0.0032781 
0.440405 - 0.0036721 
0.44® 55 - 0.0039691 
0.439865 - 0.0*2741 
0.438406 - 0.0054961 
0.436931 - 0.0062751 
0.435561 - 0.00670541 

0*5279 + 0.0005651 
0.244® 7 + 0.0®2651 
0.244506 + 0.0021061 
o. 244 ®o + 0.0030881 
0*3471 + 0.0*1961 

0*2899 + 0.0054471 
0*2316 + 0.0068441 
0*1733 + 0.0085*1 
0*1262 + 0.®00051 
0.240763 + 0.0117851 

0.239377 + 0.0222941 
0.240072 + 0.0351811 
0*3313 + 0*97681 

k 


g 21 

°31 

%L 

OOOOOOO OOOOOO 

OOOOOOOOOOOOO 

+++++++++++++ 

OOOOOOOOOOOOO 

liPillilllll 

H- H* H* H- H- h H* H- P- |S H. (I p. 

0.478588 + 0.0005371 
0.478561 + 0.0010721 
0.478498 + 0.0016061 
0.4784® + 0.0021371 
0.478295 + 0.0026661 
0.478152 + 0.0031921 
0.477988 + 0.0037111 
0.477771 + 0.0042801 
0.477573 + 0.0047431 
0.477327 + 0.0052501 
O.475616 + 0.0077381 
0.473175 + 0.0101251 
0.469953 + 0.01*4® 

OOOOOOOOOOOOO 

UtilfflU!! 

1 1 1 1 1 1 1 1 1 1 1 1 1 
OOOOOOOOOOOOO 

^mumm 

0.306807 - 0.0020941 
0.306838 - 0.0*1821 
0.306906 - 0.0062581 
O.3070® - 0.00® 211 
0.307160 - 0.01*781 
0.307341 - 0.0*081 
0.307550 - 0.*44151 
0.307843 - 0.0166061 
0.308114 - 0.®®821 

0.308456 - 0.0203281 
0.311037 - 0.0297921 
0.315002 - 0. *87591 
0.32*96 - 0. *73551 

k 

0 12 

G 22 

°32 

°42 


0.02 0.460170 + 0.0009581 0.239281 - 0.0009541 
0.04 0.460136 + 0.0019241 0.239276 - 0.0019061 
0.06 0.460062 + 0.0028981 0.239239 - 0.0028561 
0.08 0.459969 + 0.0038821 0.239191 - 0.0058061 
0.10 0.459857 + 0.0048701 0.239152 - 0.0047541 
0.12 0.459720 + 0.0058691 0.239058 - 0.0057001 
0.l4 0.459570 + 0.0068711 0.238974 - 0.0066431 
0.16 0.459372 + 0.0079811 0.238861 - 0.0076791 
0.18 0.459211 + 0.0088961 0.258757 - 0.0085281 
0.20 0.459017 + 0.0099151 0.238632 - 0.0094641 
0.30 0.457749 + 0.0150861 0.237839 - 0.0141571 
0.4o 0.456166 + 0.0203231 0.236761 - 0.0188151 
0.50 0.454244- + 0.0256231 0.235398 - 0.0234591 


-0.259283 - 0.0009591 -0.460168 + 0.00095U 
-0.239279 - 0-.0019111 -0.460133 + O.OOI8931 
-0.239245 - 0.0028701 -0.460056 + 0.0028281 
-0.239191 - 0.00383U -0.459969 + 0.0037541 
-0.239133 - 0.0047921 -0.459856 + 0.0046751 
-0.239051 - 0.0057561 -0.459727 + 0.0055871 
-0.238967 - 0.0067201 -0.459577 + 0.0064921 
-0.238857 - 0.0077821 -0.459376 + 0.0074801 
-0.238763 - 0.0086501 -0.459205 + 0.0082821 
-0.238047 - 0.0096191 -0.459002 + 0.0091691 
-0.237842 - 0.0144701 -0.457746 + 0.0135421 
-0.2367® - 0.0193171 -0.456144 + 0.0178091 
-0.235440 - 0.0*1591 -0.454202 + 0.0219951 






























TABLE I.- VALUES OF ^ ■ Concluded 


(o) M - 1.0 


k 

°10 

°20 

°30 

a 4o 

0.10 

0.12 

0.i6 

0.198675 

0.242825 

0.286975 

2.220455 + 0.0244221 
2.219995 + 0.0296811 
2.218702 + 0.0592121 
2.217279 + 0.0497461 
2.215455 + 0.0593951 
2.212373 + 0.0705181 

0.938021 + 0.0481281 
0.937372 + 0.0578241 
0.936469 + 0.0765231 
0.935663 + 0.0964131 
0.948201 + 0.1030491 
0.930797 + 0.1360241 

0.684779 + 0.&579301 
0.684637 + 0. 0700041 
0.683875 + 0.0924701 
0.677372 + 0.1170371 
0.674169 + 0.1494191 
0.679618 + 0.1661581 

0.156747 - 0.13C4801 
0.157996 - 0.1577091 
0.160954 - 0.2082051 
0.169686 - 0.2631961 
0.162175 - 0. 3U8611 
0.177212 - 0.3726801 

k 

Q 11 

°21 

°31 

°4l 

0.10 

0.12 

0.16 

O.198675 

0.242825 

O.286975 

0.277575 + 0.0053001 
0.277517 + 0.0063861 
0.277390 + 0.0083961 
0.277234 + 0. oio 424 i 
0.277012 + 0.0127341 
0.276750 + 0.0150341 

0.394836 + 0.0091531 
0.394828 + 0.0109331 
0.394583 + 0. 0146741 
0.394336 + 0.0182121 
0.393992 + 0.0222391 
0.393575 + 0.0262661 

0.480561 + 0.0126981 
0.480658 + 0.0154351 
0.480234 + 0.0205951 
0.479269 + 0.0254851 
0.479517 + 0.0312211 
0.479018 + 0.0368701 

0.180361 - 0.0271511 
0.180330 - 0.0327521 
0.181126 - 0.0436651 
0.182494 - 0.054121.1 
0.182812 - 0.0661941 
0.183990 - 0.0781701 

k 

°12 

°22 

032 

°42 

0.10 

0.12 

0.16 

0,196675 

0.242825 

0.286975 

0.052035 + 0.0013151 
0.052030 + 0.0015711 
0.051996 + 0.0020911 
0.051964 + 0.0023951 
0.051914 + 0.0031701 
0.051855 + 0.0037421 

O.I78H8 + 0.0041181 
0.178064 + 0.0049251 
0.177993 + 0.0066531 
0.177877 + 0.0082331 
0.177719 + 0.0100561 
0.177528 + 0.011878! 

0.342280 + 0.0082981 
0.342163 + 0.0099931 
0.342012 + 0.0134051 
0.34l668 + 0.0166221 
0. 341484 + 0.0203221 
0.341115 + 0.0239991 

0.227567 - 0.0137311 
0.227743 - 0.0164891 
0.227999 - 0.0221291 
0.228491 - 0.0274501 

0.228883 - 0.0335481 
0.229502 - 0.0596191 


o\ 

-J 


469c us vovn 



CT\ 

CD 


I 


VJS 2 U,- TAKES or 




f* jess aura mb o w mnxnx a 


[U - 0.?) k * 0.22] 


V 

H 

V 

."vl^ 




- 9* 



a va 

1 

0 

0.1 

0.8 

S3 

0 °:i 

?:i 

0.9 

i.ta* + o.ooaAi 
i.ran + o.ccaw 
1,68703 + 0.CD3ZQ1 
1.64233 + o,oo*4i 
1.57609 ♦ o.ooaofii 

149116 + 0. 001*1 

1.37747 + o.ooa&w 
1,225* ♦ 0.001601 
1.C5310 + 0.001351 
O.T3053 + O.OOOJ61 

0.30752 + o.owm 
000596 + 0^122301 
0.50131 + 0.0*961 
0 -M3J6 + 0.0*361 
0 J81B5 + o.oaw4i 
0.26652 + 0,01*11 

004602 + 0.017931 

o-ayrfi + 0,016011 
046451 + 0.013451 
045404 + 0.C097T1 

0 47744 + 0.011051 
043515 + 0.010961 
044630 + 0.01C3U 
045677 + 0.010551 

J + 0!^^ 

009959 + 0.0*811 

0 +01 

0.01713 + 0.000021 
0.0674a + 0.000091 
0.14765 + 0.000191 
0.23*9 + 0.000531 
0.77279 + 0.000461 

049369 + 0.000651 

0.6025a + 0.000761 
0.66118 + 0.000661 
0.60795 + 0.000791 

0 +01 
O.6C906 + 0,000221 
0.01305 + O.OOCOBi 
O.CaGX + 0.CO13CU 

0,04310 + O.OOM91 
0X6036 + 0.004831 
0,00637 + O.COftM 
0.10761 + O.OOTWl 

0.11809 + 0.000611 

0.10657 + 0.007911 

0 +01 
0,0*33 + o.ooam 

0. £30.795 + 0.0CCA31 
0,03927 + O.OOCKId 
O.O6T06 + O.Od&ii 
0.09904 + 0.002591 
O.mft + 0.005101 

0.16007 + o.ocscdi 

0.17W + O.OCWfrl 
0.16151 + 0.003901 

0 +61 

o.oocht + «■ 

0,00370 + 01 
0X0330 + 0.000021 
OXWO + 0.000031 
0.05800 + 0.000121 
0 J,T0?e + 0.000251 
+ o.oocsai 

OA«3l6 +. O.OOC031 
+ 0.000641 

0 +01 
0.000© + 01 
0.00*5 + O.OOOCAl 
O.OQ250 + o.ooarri 

0X0722 + 0X00351 

o.cojSft + o.oaizLi 
0x5169 + 0.002321 
0x377? + o.owi 
0x7330 + 0.003511 
0,00794 + q.ooftli 

0 +01 
0.00003 + 01 
0x0072 + 0.000021 
O.OQ355 + o.coo»i 
0.01073 + o.ooceSi 

O.O047< + 0.000501 

o.cii'Aj + o.ooaiAi 

0.07W3 + 0.000091 

o.n£4e + 0.002711 

0.15002 + 0,005161 

2 

0 

0.1 

o.t 

°d 

S:2 

S:J 

0.9 

0.73777 - 0.013* 
0.73367 - 0.011*1 
0.71*7 - 0.011071 
o.£?4i " 

0. 30.il - o]ot^1 

o 4 tSto + 0.003521 
04*51 + 0.005801 
040707 + O.OO3711 
0454* + 0.003361 

0.43090 + 0.003541 
041203 + 0.003*1 
0 J6136 + 0.0*661 
0.34*3 ♦ 04XAl6l 
o.aOfioa + 0.0(3501 
O.BO779 + o.ooa*i 

0.236* - 0.010351 
0.25685 - 0.010301 
0.23323 - 0.000531 
0.22706 - 0.000061 
0.20617 - 0,009671 
0,20613 - o.oogjl 

0JO3T6 - (>5*1* 

0 + Oi 

0.00734 - o.ooam 
0.02690 - 0.0C0441 
0,0ffiS3l - 0.000971 
o.iofly - o,ocn£6i 
0.15954 - O.OOrtJl 
0,11136 - 0.005251 
0,35805 - 0.00951 

0,13515 - 0X04?4i 
0,26034 - 0,005991 

0 +01 

0,0*74 + 0.0x061 
0,01668 + 0.000231 
0.040® + 0.000501 

O.O699O + 0.000631 
0.10321 + 0.001261 
0 J3T19 + 0.0*661 
0.16681 + 0.002*1 
OJJ33® + 0.002*1 
01S91 + 0.002061 

0 +01 
0.00237 - 0.000101 
0X0««3 - O.OOftli 
o.oaoMi - 0.0009U 
0.05491 - 0.001551 
O.CSlft - 0.002201 
O.CSB33 - 0.0030^1 
0.003*) - 0. 005^91 
O.OgAO - O.OCWSi 
O.OWOJ - 0.005131 

0 +01 
0.00007 + 01 
O.OQllb - O.COOOH 
0.00370 - O.COC091 
0.01730 - O.OOCC71 
0.TOS1 - 0.0006U 

0,07643 - O.COU.71 
002643 - 0.000941 
oafiua - 0x02781 
0.0006 - 0X0531 

0 + * 
0.00001+ 01 
o.ooow + o^oocai 

0.00568 + 0.000*1 
0.01118 + 0.00*41 
0.02380 + 0.000311 
0.**2 + 0.0006* 
0.06174 + 0.0*0* 
0.11713 + 0.0*451 
0.13633 + 04x0671 

0 +01 
0.00008 + 01 
O.OCO57 - O.COOOfei 
o.ooaA - 0X00001 
0.0053? - 0.000851 

0.0CL2EU - O.OOOJti 
0.0*60 - 0.001091 
0.04002 - 0. 001011 
0.03030 - OXQ2391 

0X6000 - 0.003021 

5 

0 

0.1 

0.8 

°o* 

0.5 

0.6 

°o.l 

0.9 

0,14132 - 0.0*841 
043531 - 0.0*811 
0.45164 - 0.0*741 
o. 4 sob - o.oAsei 
0405* - 0.0*431 
0.36130 - 0.0*191 
0.35*2 . 0.0035ri 

0.51439 - o.otgiti 
0.26431 - 0.002901 
0. 19102 . 0.002m 

04T702 - 0.003631 
047463 - 0.0038U 
0.40738 - 0.007721 
045503 - 0.005771 
0.437*0 - 0.C05351 
041311 - 0.005051 
o.53i£o - 0.0*671 
0.54066 - 0.0*171 
0.286* . 0.005301 
0.20795 - o.ooejn 

-0.25755 - 0,010731 
X. 256ft - 0.010611 
-0.85275 - 0.000311 
-0.aa539 - 0.000831 
-0.21770 - 0.009C31 
-0.80m - 0.009691 
-0.1900a - o.oonai 
-0.16965 - 0,007661 
-o.Uaja - o.ocwi 

-O.lC^. - 0.0CA6&1 

0 +01 
0.00438 - 0.000051 

0.01T26 - 0.00*91 
0.03T8S - 0.00*21 

0.06460 - 0.01x711 
0.09337 - 0.001631 

002697 - 0.001391 
005415 - 0.001701 

006916 - 0. 0*661 
0.13554 - 0.001711 

0 +01 
0.0*75 - O.OCOO61 
O.CQfljo - 0.00*51 
0.*C+> - 0.000301 
0.06993 - 0.000061 
0.1D3S8 - 0.0*261 
0.13758 - 0.0016B1 
0 JL6^2 - 0.002*1 
0.18317 - 0.002*1 
0.16*2 - 0.002061 

0 +01 
-0.00 256 - 0.00Q1U 
-0.00951 - o.oocto 
-OX«059 - O.COOSd 
-0.05M6 - 0.001571 
-0.05143 - O.COSfil 
-0,06^1 - 0.000091 
-0.00512 - o.ami 
-o.op2i - o.octoi 
-o.o050r - 0.CX9791 

0 + w 

0.0X04 + * 
0J»069 - 0.000011 
0.03340 - 0.000*1 

0.01034 - o.ooorn 

0.02364 - 0.000661 
0.*367 - 0.0005* 
0.07533 - o.ooofei 
0.10&6 - 0.001191 
0.12399 - 0.0*391 

0 + * 
0.00003 + 01 

0.00013 - O.OCOOll 
0.00369 - 0.000*1 
0.01119 - 0.00*41 
0.0238 b - 0.000511 
o.*94* - 0.0006* 
;o,ceiT9 - 0.0*0* 
0.11725 - 0.0*451 
0.13*2 - 0.0*671 

0 +01 
xxoooa + oi 

-0,00057 - 0.000021 
-0.001ft - 0.000081 

-0.00397 - 0.000251 
-0X1206 - 0x00581 
XXA63 - 0,0am 1 
-0x4073 - O.OOlftl 
X. 05057 - 0.002ft i 

-0x6795 - 0,003071 

1 



0.30056 - 0.02*01 

0.50W1 - 0.023291 
0.30531 - 0.021991 

O.fiOTO - 0.021971 
0.26572 - 0,020531 
0x6009 - 0.01*01 
0.AT65 - 0.ca7»l 

0.C210T - 0.016001 
0087* - O.QUttl 

0.1*95 - o.ooyrn 

mm 

pi! ! i'$k 

[ 322 ! : S2E 

r 7?*?. * : rriit-T 

e sa • 

0 +01 
o.ocfti + 0.0001H1 

O.OOft0 + O.COQfti 
0.0607* + 0.001191 
0.05ft7 + 0.002051 
0.05256 + 0.005011 
0 A 6966 + O.OftOQl 

oxete + 0.00W61 

0,09208 + 0.005341 
0,00539 + 0.004911 

1 *iLjC * 

: j ^PS 

0 + * 

-0.0*56 + 0.000101 
-oxaT® + o.cootoi 
-0,05952 + o.coMBi 
-0.06716 + 0.0*5* 
-6.09915 + 0.002221 

-0.13189 + 0.002951 
-0-16323 + 0.003581 

-0473*3 + O.OCBB1 
-OJ6169 + O.OC0611 

0 +01 
0.00002 + Od 
0,00030 + 0X00021 

o.ocnfrr + o.ooam 
O.OO36B + 0.000001 
o.cnjo? + 0x00731 

0,00300 + 0.0QW41 
0.C4147 + 0X02501 
0.055*4 + 0.0Q3421 
0.06917 + 0,0(25961 

0 +01 
0.00005 + Oi 
0.00<*9 - 0XCO041 
0.00639 - O.OOCa7± 
o.wpfi - 0.000321 
O.ca ten - o.ooini 
• 0.05809 - 0.002321 

0.03500 - 0.003641 
0.07606 - 0.003301 
a.o«3t - o.oofiku 

s N S+ 

IB 

f. m J pi 


NACA TN 3694 





J 


a hi in.- fuczs or re eootuss tactko 


(k - o.ji'it - o.s^ 


I Control point 1 I Control point 2 I Oarteol point 3 Ooctrol joint 4 Control point 3 I Control point 6 I Control point 7 I Co nt r o l point 8 | Control joint 9 


1 0 -0.21173+ 0. 

0.1 - 1 . 32029 + 0 . 

o.t 4.07351-0. 

0.3 -1.26418+ 0. 
0.3 -0.22782+ 0. 
0-3 - 0 . 00 . 604 . 0 . 

0.6 -0, 1*939+6. 
0.7 -0.03073+0. 
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Figure 5*- Variation of moment with reduced frequency k for a delta wing 
oscillating in pitch about its midchord. A - k; M = 0. 
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variation or nru ana morasnx warn reaucea irequency k ror a recranguiar wing oscii- 
U in pitch about its midchord. Comparison of results of multiple -line method with results 
[rface-loading method. A ° 2; M «* 0. 






